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Abstract 

A biunimodular vector of a unitary matrix A £ U (n) is a vector v £ T n C C n such that 
Av £ T n as well. Over the last 30 years, the sets of biunimodular vectors for Fourier matrices 
have been the object of extensive research in various areas of mathematics and applied 
sciences. Here, we broaden this basic harmonic analysis perspective and extend the search for 
biunimodular vectors to arbitrary unitary matrices. This search can be motivated in various 
ways. The main motivation is provided by the fact, that the existence of biunimodular 
vectors for an arbitrary unitary matrix allows for a natural understanding of the structure 
of all unitary matrices. 


Note 

As we have learned from the referees, the existence of biunimodular vectors for an arbitrary 
unitary matrix has been recently proved in this journal by Idel and Wolf (see [36| ) basing on 
the results of [9] in symplectic geometry due to Biran, Entov and Polterovich. In this paper 
we provide an independent account of this existence phenomenon and push a bit further the 
theory of biunimodular vectors. 


Outline 

In Section 1 we describe the background for this research. The next two sections discuss 
the structure of all unitary matrices via decompositions based on biunimodular vectors. In 
Section 4 we state the problem of finding biunimodular vectors in several equivalent terms. 
One of them is a maximal torus decomposition of SU(n ), that can be viewed as a close relative 
of the spectral theorem. In the following sections we justify that all matrices in U( 2) and U (3) 
possess biunimodular vectors and provide two closed formulas for U (3). One of them being an 
analogue of the standard Euler-Tait-Bryan decomposition of 0(3). In Section 8 we provide a 
simplified novel method to synthesise all matrices in [7(2"), that enjoys an elementary proof 
for all n £ N. Section 9 is devoted to the study of biunimodular vectors in the general case 
from the Lie groups perspective. Using basic tools we show that in arbitrary dimension n, the 
set of matrices possessing biunimodular vectors contains a nonempty open subset of U(n). 
We also obtain that, up to normalization, the set of biunimodular vectors is finite for almost 
every unitary matrix. Finally, we employ a simple alternating projection algorithm that for 
a numerically given unitary matrix produces vectors, that are biunimodular within a set 
precision. We test this algorithm to provide numerical evidence indicating its effectiveness 
in finding such near biunimodular vectors for random unitary matrices with dimensions up 
to 100. We also apply the algorithm to Fourier matrices, and compare the search results to 
the known solutions. 
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1 Background 


The paper explores the notion of a biunimodular vector that traces back to Gauss. The 
classic example of such vectors are Gauss sequences given for k = 0,1,..., n — 1 as 


u k = ce^ (Afe2+/ifc) 


n odd, 


u k = ce 


-(^fc 2 +/ifc) 


n even, 


where A, p £ {0, 1,. .., n — 1} with A relatively prime to n € N and c € T = {z G C : |z| = 1}. 
These vectors are unimodular, that is, [uk\ = 1 for all k = 0, — 1 and have an 

interesting property. Their Fourier transform F n u, defined by aplying to u the unitary nx n 
matrix F n with entries ( F n )j^ = ~^r l e ~ 2m "A is unimodular as well 1 

Following Haagerup ([31]) and Saffari (HU) we attribute the starting point of the general 
theory of biunimodular vectors to Per Enflo. In 1983 he has asked whether for prime n the 
only unimodular vectors of length n whose Fourier transform is also unimodular are Gauss 
sequences. While the answer is affirmative when n < 5, a computer search done by Bjorck 
for n = 7 provided a surprising counterexample 


(1; 1, 1, e ) 1) e ) e ) 


with 6 = arccos(—|) yielding e l9 = — | + i^-. These findings, published in 1985 ([TO]), 
have been generalized in HU and later the term “bi-unimodular sequence” has been coined 
by Bjorck and Saffari to describe a unimodular finite vector, whose Fourier transform is 

unimodular as well (see nan- 

The Bjorck sequences provided in hh are biunimodular vectors of a prime length p. 
For p = — 1 (mod 4) their coefficients are either 1 or e l9 with 9 = arccos For p = 1 
(mod 4) their coefficients are either 1, e lv or e~ lri with rj = arccos ^d +1 . The importance 
of Bjorck sequences has been underlined in [5], where it was shown that for these sequences 
(contrary to Gauss sequences) the discrete narrow band ambiguity function has an optimal 
bound. However, Gauss and Bjorck sequences provide only a glimpse at the structure of 
biunimodular vectors for F n . Even in the case p = 7 there is another such vector found by 
Bjorck and Froberg that is related neither to a Gauss nor to a Bjorck sequence (see m or 


Moreover, in [T5] Bjorck and Saffari proved that if n is divisible by a square, then there 
are infinitely many (continuum) biunimodular vectors with the leading entry 1, and they 
conjectured that for all square-free n the number of such biunimodular vectors is finite. 

Conjecture 1.1. (Bjorck, Saffari) If n £ N is not divisible by a square, then the set of 
biunimodular vectors (with the leading entry 1) for the Fourier matrix F n is finite. 

An easy observation is that a unimodular vector u = (uq,ui, ... ,u n — i) is biunimodular 
for F n if and only if its cyclic translations are orthogonal. This orthogonality relation is 
called zero autocorrelation and means that 

n— 1 

y; u k Uk+i = 0 for i = 1,2, ... ,n — 1, (1) 

fc=o 

where the index k + i is taken (mod?r). 

1 As a matter of fact, the Fourier transform of a Gauss sequence yields another Gauss sequence. 
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The idea of Bjorck to search for biunimodular vectors u for F n relies on setting Xk = 
where k = 0,1 ,..., n— 1 with u n = u 0 and transforming 0 to the following set of equations: 

Xq + %i + ■ ■ ■ + x n —i = 0 

X 0 X\ + X\X2 + ■ ■ ■ + X n -iX 0 = 0 


X 0 Xi ... x n -2 H - 1- X n -iX 0 ■ ■ ■ X n -3 = 0 

XqXi . X n -1 = 1 


( 2 ) 


The unimodular solutions of the above set of equations are in one-to-one correspondence to 
biunimodular vectors for F n (with the leading entry 1). The general complex solutions of 
([2]) are called cyclic n-roots and became a research area of independent interest. In the mid- 
1980s Arnborg has asked Davenport to establish whether the set of cyclic 6 -roots is finite, 
leading to the popularisation of the cyclic n-roots problem in [22] . Backelin proved in [2] that 
for n divisible by a square the number of cyclic n-roots is infinite. For 2 < n < 8 all cyclic 
n-roots were found by Bjorck and Froberg (see p~2l IITTI ITTj 1 by using Grobner basis techniques 
to solve the system |2]). Faugere has found cyclic 9 and 10 roots by improving the search for 
the Grobner basis (see m and [25]). These results have been confirmed and extended by 
polyhedral homotopy continuation methods applied for solving 0 as a benchmark problem. 
The numbers of cyclic n-roots for a square-free n between 2 and 14 is listed belowj^j 


n 

2 

3 

5 

6 

7 

10 

11 

13 

14 

15 

cyclic n-roots 

2 

6 

70 

156 

924 

34,940 

184,756 

2,704,156 

8,464,307 

? 

biuni vectors 

2 

6 

20 

48 

532 

? 

? 

53,222 

? 

? 


The table also contains information about the number of unimodular cyclic n-roots that 
yields the amount of biunimodular vectors for F n . The numbers up to n = 7 were given 
by Bjorck, Froberg and Haagerup, who inspected the obtained cyclic n-roots to check for 
unimodular solutions, establishing the real benchmark for solving the cyclic n-root problem 
(see 1 12] and [34]). The number of biunimodular vectors for n = 13 has been obtained by 
Gabidulin and Shorin in [301 . 

Basing on the number of cyclic p-roots for prime p up to 7 Froberg concluded that the 
amount of cyclic p-roots should be ( 2p ~ L 2 ) ■ This conjecture has been confirmed by Haagerup 
in [55] under the condition that the roots are counted with multiplicities. 

Theorem 1.2. (Haagerup) For prime p the number of cyclic p-roots counted with multiplic¬ 
ities is ( 2 pZi) ■ 

The paper of Haagerup contains an elegant direct argument showing that the number 
of cyclic p-roots is finite for prime p, confirming Conjecture 0 in the prime case. His 
more advanced reasoning allows to count the cyclic p-roots with multiplicities, imposing a 
question whether for all primes cyclic p-roots have multiplicity 1 and providing an additional 
motivation to verify the number of distinct cyclic 13-roots. 

A slightly different approach towards finding biunimodular vectors for F n has been taken 
by Gabidulin and Shorin. In [30] they multiplied the system (JTJ) by the product JIfc=o u k *o 
obtain a system of equations that has been treated with Grobner basis techniques as well, 
providing the mentioned number of biunimodular vectors for F 13 . 

The theory of biunimodular vectors may be hard to follow due to their popularity stem¬ 
ming from applications in several fields of signal processing. This popularity can be mea¬ 
sured by the amount of various names given to the same objects by many researchers working 

2 For 2 < n < 11 the amount of roots has been verified by various methods, for n = 13 and 14 the number of 
cyclic n-roots is claimed only by Li and Tsai (see [13). Moreover, they count only isolated roots, so the number 
of cyclic 14-roots still has to be verified as finite or not. 
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within the area. Gauss sequences are often called Zadoff, Chu or Wiener sequences. Unimod- 
ular vectors are described as constant amplitude, phase-shift keyed or polyphase sequences. 
Vectors whose Fourier transform is unimodular are named as perfect, having optimum cor¬ 
relation or zero autocorrelation. Therefore, biunimodular vectors are studied as CAZAC 
(constant amplitude zero autocorrelation) sequences, PSK (phase-shift keyed) perfect se¬ 
quences, polyphase sequences with optimum correlation and so on (see [7] by Benedetto 
et al. for a better account of this nomenclature and a list of relevant papers that can be 
extended by adding an early reference HD- 

Since full understanding of biunimodular vectors would allow to resolve the circulant 
Hadamard matrix conjecture, and is just a tip of the iceberg being a classification of all 
Hadamard matrices, the complexity of the topic is apparent. 

It is important to notice, that the research on biunimodular vectors is not confined to 
the Fourier transform F n on the cyclic group Z„. Indeed, real biunimodular vectors for the 
Fourier transform on the group {hf) 71 are investigated under the name of bent functions. Here 
n must be even and the number of bent functions is known for n < 8. The latest count is due 
to Langevin and Leander who proved the number to be 99270589265934370305785861242880 
in the case n = 8 (see [55]). 

In [51] it has been proved that biunimodular vectors exist for the Fourier transform on 
any finite abelian group. This raises the question of existence of such vectors for an arbitrary 
unitary matrix and motivates the following 

Definition 1.3. Let A be annxn unitary matrix. A vectorv £ C n is called a biunimodular 
vector for A if both v and Av are unimodular (i.e. have all entries of modulus one). 

One can argue that the above notion has surfaced already in works on the problem 
of mutually unbiased bases (MUB). Indeed, the problem can be brought down to finding 
all biunimodular vectors for a given complex n x n Hadamard matrix and checking for a 
subcollection of these vectors that forms an orthonormal basis of C n . This approach has 
been used in [33] by Grassl for Fq and by Bengtsson et al. for some other 6x6 Hadamard 
matrices (see m- In [15] Brierley and Weigert conducted an extensive search on all known 
Hadamard matrices in dimension 6 providing an impressive study of their biunimodular 
vectors leading to an approximate solution of the MUB problem in this dimension (for lower 
dimensions accessed via biunimodular vectors see IISD- 

in the following sections we explain that the existence of biunimodular vectors for an 
arbitrary unitary matrix allows to understand the structure of unitary matrices in fairly 
simple terms. 


2 Synthesis of U{n ) 

Unitary n x n matrices form the unitary group U(n) and are a basic notion in mathe¬ 
matics, physics, chemistry and engeneering. These matrices can be obtained in various 
standard ways: Gram-Schmidt process, Householder reflections, Hurwitz parametrization, 
Cayley transform and via Hermitian matrices using the matrix exponential. In 1952 Mur- 
naghan described a convenient parametrization of the unitary group (see [42j){^| This has 
been followed by the work of Reck et al. who in 1994 proved that any unitary matrix can be 
obtained as a sequence of consecutive beam splitter transformations, that can be conducted 
experimentally in the laboratory using optical devices (@5j). In the current century the 
research on finding the ways to generate unitary matrices has intensified. Nemoto described 
the special unitary group SU (n) in |43j basing on an earlier result [46] of Rowe et al. Tilma 
and Sudarshan provided an Euler angle parametrization for SU(n) and U(n) (see [501 l5l] ). 
Meanwhile Di^a in [25] has shown another description of U{n). A recursive parametrization 


3 Another such parametrization has been provided in 1982 by Di1,a (see [24]'). 
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of unitary matrices has been given by Jarlskog in [35] and a composite parameterization was 
done in [35] (see mm and @ 8 ] for more details on these developments). 

The search for biunimodular vectors for an arbitrary unitary matrix can be motivated by 
the following recipe to generate unitary matrices. Consider n £ N and a set of n 2 complex 
parameters {djk £ T : j,k = 1,2,, n}, where T = {z £ C : \z\ = 1}. In order to synthesize 
a unitary matrix A n £ U(n) from these parameters, it is convenient to organize them as in 
the table below 


0\n 
0-2 n 
03n 


Onl O n 2 O n 3 . . . O nn 


an 

012 

Ol 3 

023 

a 21 022 

031 032 033 


and run the following recursive procedure using the Fourier transform (Fi)j^ = e 
and its conjugate transpose (. F*)j^ = ^e 2l, r . 

The initial unitary matrix Ai is given as [an]. And, in general, for 2 < l < n 



ail 

0 . 

. 0 


1 

0 ... 

0 



'1 

0 . 

0 

Ai = 

0 

ai2 ■ 

. 0 

Ft 

0 

Ai-i 



Ft 

0 

an ■ 

0 


_ 0 

0 . 

au_ 


0 





_0 

0 . 

ai~i i_ 


with an obvious block matrix notation. 

The above procedure is not one-to-one, meaning that different sets of parameters may 
give the same matrix. As to the question if it is onto, the simplicity of the above scheme 
can raise some doubts whether it can reveal the structure of an arbitrary unitary matrix. 
Nevertheless, in the next section we shall prove that all unitary matrices can be obtained in 
this way, provided that all of them possess biunimodular vectors. 


3 Analysis of U(n) 

We have the following basic 

Proposition 3.1. A unitary matrix A £ U{n ) has a biunimodular vector v if and only if 



Wl 

0 . 

. 0 


1 

0 

... 0 



Vi 

0 . 

. 0 

A = 

0 

w 2 ■ 

. 0 

F n 

0 


B 


K 

0 

V2 ■ 

. 0 


_ 0 

0 . 

• w n _ 


0 





_0 

0 . 

• ty 


where v = (t>i, v 2 , ■ ■ ■, v n ) £ T", Av = w = (wi,w 2 , ■ • •, w n ) £ T n , B is some unitary matrix 
in U(n — 1) and F n is the Fourier transform with the adjoint F*. 

Proof. We shall use the following notation. For an arbitrary vector u £ C n by D u we 
shall denote the diagonal matrix with u on the diagonal. We also consider the stabilizer group 
of u, defined as Fix(u) = {S £ U{n) : Su = u} and let 1 = (1 ,1,..., 1), e = (1, 0,..., 0) 
both in C". 

To conduct the proof we observe, that there exist unimodular vectors v and w such that 
Av = w iff AD V 1 = D w 1, what is equivalent to saying that D^ADy £ Fix(l) with v, w £ T". 
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The description of the group Fix(l) is related to the group Fix(e). Indeed, since F n e = 
^=1, it is easy to see that M £ Fix(l) iff F*MF n £ Fix(e). 

Therefore, we see that A has a biunimodular vector v iff F*D^AD v F n £ Fix(e). Due to 
the obvious description of Fix(e) the proposition follows. □ 

Remark 3.2. Clearly, if we can find a biunimodular vector v for A, then for every c £ T the 
vector cv is also biunimodular for A. Therefore, we can concentrate on biunimodular vectors 
with the leading entry v\ equal to 1. This explains why Proposition \3.1\ proves, that every 
unitary matrix can be obtained via formula 0 iff every unitary matrix has a biunimodular 
vector. 


Remark 3.3. The Fourier transform F n enters the above scheme as a generic unitary matrix 
with the property that F n e = -^1. We need to clarify that, in the above proposition, F n 
can be replaced by any unitary matrix F such that Fe = -^= 1. However, clearly the easiest 
description of such F is given by F = F n M with M £ Fix(e). 

In Section [To] we shall present an algorithm that for an explicitly given unitary matrix 
allows to search for its biunimodular vectors within a set precision. Therefore, in practice, a 
numerically given unitary matrix can be analysed to the point when all parameters {ojk £ 
T : j, k = 1, 2,..., n} are found, that allow to synthesize the matrix back via formula (13]). 

On the other hand, as pointed out by the referees, one has the following crucial 

Theorem 3.4. (Idel, Wolf + Biran, Entov, Polterovich) Every unitary matrix has a biuni¬ 
modular vector. 


This settles that scheme |3]) for building unitary matrices gives all of them. The scheme 
has been discovered independently by De Vos and De Baerdemacker without the notion 
of biunimodular vectors. In |23j they conjectured that every unitary matrix U enjoys a 
LAR decomposition: U = LAR , where L and R are unitary diagonal matrices and A is a 
unitary doubly stochastic matrix, meaning that the row sums (and thus columns sums) of 
A are all equal to 1. This, clearly, is equivalent to saying that every unitary matrix U has 
a biunimodular vector and has been proved by Idel and Wolf in [36| . The proof relies on 
showing that the action of U on the torus T n defines a Hamiltonian symplectomorphism of 
the corresponding Clifford torus embeded in the complex projective space. The rest follows 
from 9j, where Biran, Entov and Polterovich used Calabi quasimorphisms to prove that the 
Clifford torus cannot be displaced from itself by a Hamiltonian symplectomorphism. Idel 
and Wolf poin t also to the work [2D] of Cho on Floer cohomology that one can combine with 
Theorem 3.4 to conclude, that every n x n unitary matrix has at least 2 Tl ~ 1 biunimodular 
vectors. 

In the following section we shall examine the existence of biunimodular vectors from a 
few different angles. 


4 Biunimodular problem 

The biunimodular problem is to find (or prove the existence of) a biunimodular vector for 
a given unitary matrix. With the notation explained below we can restate this task in 
equivalent terms. 

Theorem 4.1. Let A £ U(n). The matrix A has a biunimodular vector if and only if either 
of the following holds 

(a) A £ V ■ Fix(l) • T>\. 

(a’) A! £ V ■ E(l) • V. 

(b) ||A||oo—).! = n. 
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(c) sup„ eT „ |tt a (w)| = 1 . 

(d) H(T n ) n T" / 0. 

Proof. 

(a) Here, V C U(n) denotes the subgroup of diagonal matrices with diagonal entries 
in T, and D\ C T> is the subgroup of matrices with the first diagonal entry equal to 1. 
Fix(l) = {S £ U(ri) : S 1 = 1}, where 1 = (1,1,..., 1) G T n . A £ V ■ Fix(l) • V\ means that 
A can be written as a product of three matrices from the indicated subgroups of U ( n). By 
Proposition |3.1| and Remark |3.2| this is equivalent to the existence of a biunimodular vector 
for A. 

(a’) Here, A! = AH with A ra = det(H), so A' £ SU(n). V = T>r\SU(n) and E(l) denotes 
the subgroup of SU(n ) that has 1 as an eigenvector. Clearly, (a’) is a mere translation of 
(a) into the special unitary setting. 

(b) Here, || A^^ = sup || u || Qo=1 PuHr is the norm of A : (C n , || ■ IU) -» (C n , || • ||i), 

where the usual ^ p -norm of u £ C n is given by ||w|| p = (X^k=i \ u k\ p )* and ||u||oo = sup fe \uk\- 
For every u £ C n we have ||«||i < -y/n||u|| 2 < ro||zt||oo- Therefore, for every A £ U(n ) we 
have 

||Hu||i < y/n\\Au || 2 = Vn\\ u h < «IMIoo> ( 4 ) 

proving that ||^4||< n. Thus, if A has a biunimodular vector v, then ||Hu||i = nH'uHoo 
showing that HHHoo^i = n. 

Conversely, if IjHHoo^i = n, compactness of the f 00 -ball implies the existence of a vector 
u £ C n with ||u||oo = 1 and ||Hzt||i = nHuHoo- As the estimate 0 shows, this equality holds 
if and only if 1111! = and ||u || 2 = \/^IMIoo! meaning that both u and Au must 

be a constant multiple of a unimodular vector. Since ||rt||c» = 1 and ||Hu||i = n, both u and 
Au must be unimodular. 

(c) Here, tta(u) = n t{Au), where tt(u) = JIfe=i u k for u £ C n . 

Let A £ U(n), u £ T n and Au = w. We have 



so |7 Ta(w)| < 1. Thus, if A has a biunimodular vector v, then |7 ta('c)| = 1 showing that 
su Pii£ T " l 7r ^(w)| = 1. 

Conversely, if sup ueTrl \tta(u)\ = 1, compactness of T” implies the existence of a vector 
u £ T n with |7Ta(u)| = 1. In light of (|HJ) this means that 1 = IIfc=i l w fc | 2 = de. 

w = Au must be unimodular. 

(d) Clearly, a vector v is biunimodular for A £ U(n) iff Av £ H(T n ) D T n , so the 
equivalence in this case is trivial. □ 


Remark 4.2. Part (a) of Theorem 4-1 is equivalent to saying that A £ 2? • Fix(l) ■ 2?. Thus, 
Theorem 3.4 is equivalent to the decomposition U(n) = T> ■ Fix(l) • T>. Unfortunately, there 
is no general theory explaining when two subgroups H, H' of a Lie group G yield the decom¬ 
position G = H ■ H' ■ H. Moreover, the Cartan decomposition and Kostant decomposition of 
any semisimple Lie group G, that both are given in this form, indicate potential difficulties 
in developing such a general explanation. 

For part (a’) included in the above characterization, we would like to notice that V is 
a maximal torus subgroup of a compact simple Lie group SU(n). Thus, Theorem 3.4 is 
equivalent to the following maximal torus decomposition: 


SU{n) = V ■ E(l) • V, 


with dim(SU(n)) = dim(E(l)) + 2dim(iy). This allows to conjecture, or ask, whether every 
compact simple Lie group G admits a maximal torus decomposition G = T ■ H ■ T, with a 
maximal torus T of G and a subgroup H C G such that dimG = dimiJ + 2dimT. 
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Regarding part (b), as in 131}/ . an interpolation argument can be used to conclude that 
||^4-||oo—»-i = n is in turn equivalent to saying that ||A|| p _» g = n^~p in the region 1 < q < 2 < 
p < oo. This further amplifies the interest in biunimodular vectors, that are precisely those 
vectors, where the norm ||A|| p ^ g = is attained, in the indicated range of p and q. 

Part (c) of the theorem is equivalent to sup ueD „ |7 t,a(u)| = 1, where B = {zeC:|,z|<l}. 
Thus, the biunimodular problem can be viewed as an extremum problem for the modulus of a 
polynomial : D n —> D. 

Part (d) of the characterization provides a deceptively simple geometric interpretation of 
the problem: Describe the intersection of two tori. This geometric angle has been exploited 


to prove Theorem 3-4 


With the above characterization we close the motivational part of the paper and present 
our findings on the biunimodular problem. 


5 U (2) 

Obtaining biunimodular vectors of a given unitary 2x2 matrix is an easy task. 
Proposition 5.1. Every matrix in U( 2) has a biunimodular vector. 

Proof. Every matrix A £ U{ 2) can be written as 


A = 


x 

-zy 


V 

zx 


with z £ T and x,y £ C such that \x\ 2 + \y\ 2 = 1. Finding a biunimodular vector (l,d) 
for A amounts to solving a system of two equations: \x + yd\ = 1 and \y — xd\ = 1. Since 
\x\ 2 + \V? = \d\ = 1 the system is equivalent to the equation R e(xyd) = 0. Therefore, we see 
that if xy = 0 then every vector (1, d) with d £ T is biunimodular for A. While for xy ^ 0 
the matrix A has precisely two biunimodular vectors with the leading 1 entry: 


(1 


(1 


( 6 ) 


where a = arg x and /3 = arg y. □ 

As we indicated in Section [2j the existence of biunimodular vectors for every matrix in 
17(2) allows for writing an alternative formula describing 17(2). 

Corollary 5.2. 


j 


a(l + z) ac(l — z) 
6(1 — z) 6 c(l + z) 


: a, b, c, z £ T 


(7) 


Proof. Let A £ 17(2). Since A has a biunimodular vector, by Theorem 4.1 we have that 
A £ V ■ Fix(l) • T> i. Clearly, in this case 


Fix(l) = 
and therefore 


ri 

'1 

1 ' 


'1 O ' 


'1 

1 ' 

nrl f 1 [ 

12 

1 

-1 


0 2 


1 

-1 



1 + z 1 — z 
1 - z 1 + z 


: z £ T 




a 0 


0 b 



1+z 1-z 

l-z 1+z 


with a,b,c,z£ T as desired. 


□ 

To close this section we point out that the above observations follow, as well, from 


writing an arbitrary matrix in 17(2) as 
ip £ [ 0 , 27t] and noting that 


cos p 
- sin p 


a 0 

0 6 

siny> 
cos p 


cos p 
— sin</? 


smyj 

cosv? 



'1 O' 


0 c 


with some a,b,c £ T, 


_ e~ itp 

'1 o' 


2 

0 -i 



1+z 1—2 

1-2 1+2 



'1 O' 


0 i 


with 


2 = e i2ip 
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6 U (3) 


Some elementary calculations allow to conclude that all orthogonal matrices in 0(3) have a 
biunimodular vector. Another easy case is given in the following. 

Example 6.1. For A £ U( 3) with at least one zero entry we shall explain how to exhibit 
its biunimodular vectors. By interchanging rows and columns of A we can assume that there 
is a zero entry in the top right corner of A. Moreover, we can multiply rows and columns 
of A by unimodidar constants in such a way that the outcome shall have only real entries. 
Therefore, we can assume that A is given as 

cos a sin a 0 
a; sin a — 2 cos a y 
y sin a —y cos a — x 

with a £ [0, 2n\ and x,y £ R. such that x 2 + y 2 = 1. A simple check reveals that such a 
matrix has at least four distinct biunimodular vectors with the leading 1 entry: 


(1 




(1, i, -e ia ), 


(1 e" ia ). 


Moreover, if A has only one zero entry, then there are no other biunimodular vectors for A. 
It is also clear, that when A has at least two zero entries, then at least one of its entries has 
modulus 1 leading to a trivialization of A and a continuum of biunimodular vectors. 

Remark 6.2. Whenever we talk about the number of biunimodular vectors for a given unitary 
matrix, we concentrate only on vectors with the leading entry equal to 1, as it was indicated 
in Remark \3.2\ The above example explains the structure and the number of biunimodular 
vectors for A £ U( 3) in the case when A has at least one zero entry. In the case when all 
entries of A £ 17(3) are nonzero, we have only found examples of matrices with the number 
of biunimodular vectors equal to 4,5 or 6. 

In order to exhibit the method allowing to count biunimodular vectors of a given matrix 
A £ 17(3), let us consider a vector u = u xy = (l,e“,e ly ) and the square Q = [— 7 r, 7 r] 2 . The 
idea is to look at three regions Rj = {(a :,y) £ Q : \Au{j)\ > 1 }, j = 1,2,3 and check the 
intersection of their boundaries. 

For the Fourier case A = F 3 the corresponding three regions are given in FigjTJ The points 
where the boundaries of these regions intersect correspond to biunimodular vectors of A. 
From the unitarity of A it follows that every intersection point of two such boundaries belongs 
to the boundary of the third region as well, allowing for a visual count of biunimodular 
vectors. The next figure shows these regions for three exemplary unitary matrices providing 
the mentioned count of biunimodular vectors (conducted and plotted with Mathematica). 

In general, since A is unitary these closed regions cover Q, and therefore, any of these 
regions must intersect with at least one of the two other regions. Unfortunately, it is hard 
to argue that at least one pair of the regions intersects in such a way that their boundaries 
also intersect, closing a natural way to prove Theorem |6.3[ This obstacle boils down to the 
lack of proof for a classification of the regions Rj, that essentially can be simplified to 

E st = {{x,y) £ Q : scos(x) + tcos(y) + cos(x — y) > 0}, s,t > 0 
and are hard to study. 

A working proof that all A £ U (3) have biunimodular vectors, that we give below, relies 
on findings of Section |HJ that are included in Appendix A. 

Theorem 6.3. Every matrix in 17(3) has a biunimodular vector. 





Figure 1: The regions Rj, j = 1,2,3 for the Fourier case F 3 . 





Figure 2: Counting biunimodular vectors in dimension n = 3. 
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Proof. In Appendix A (Observation |A.3| we prove that 
U{ 3) = {diag(Ai, A 2 , A 3 )r (a>/ 3 > 7 ; 2 ) diag(l, A 4 , A 5 ) : a, 0, 7 € [0, 2?r]; z, Xj £ T, j = 1,..., 5}, 
where 


,z) ~ 


cos a — sin a sin 7 sin a cos 7 

— sin a sin 0 z cos 0 cos 7 — cos a sin 0 sin 7 z cos 0 sin 7 + cos a sin 0 cos 7 
sin a cos 0 z sin /3 cos 7 + cos a cos 0 sin 7 z sin 0 sin 7 — cos a cos 0 cos 7 


Thus, in order to prove that every matrix in U{ 3) has a biunimodular vector, it is enough to 
show it for the matrix 

Therefore, we want to find x,y £ R such that T/ a ,p^ >z ) applied to the vector (l,e lx ,e ly ) 
yields a unimodular outcome v £ T 3 . The first entry of this outcome is given by cos a + 

cos gl sin q 

w sin a with w = —e lx sin 7 + e ly cos 7 . By applying . to (1, w), we see that 

sin ol cos cy. 

| cosa+u>sina| = 1 iff |u>| = | sin a—wcosa| iff sin a—w cos a = e lm w for some m £ R. Thus, 
we can add to as another variable to our problem and require that sin a — w cos a = e lm w. 

Then, the second entry of v is given as V 2 = uz cos 0 — e lm w sin 0 , where u = e lx cos 7 + 
e ly sin 7 . And V 3 = uz sin/3 + e lrn wcos0. In short, 


V 2 


. V 3 . 



cos 0 — sin 0 

sin 0 cos 0 


0 

Am 


cos 7 sin 7 
— sin 7 cos 7 


o*y 






( 8 ) 


and we need to find x,y 7 m £ 
with w = —e lx sin 7 + e ly cos 7 . 


R such that |u 2 | = |u 3 | = 1 and sin a — wcosa = e lm w , 
Since Q reads U m (e lx , e ly ) S T 2 with U m £ 1/(2), we can 
recognize the task of finding x , y,m as a biunimodular problem with a certain twist. 

By looking at © we can simplify the above problem. Indeed, by Proposition |5.1[ for 
every to £ R we can find at least two vectors (1 ,±ie lv>m ) such that U m (l, ±ie xy,m ) £ T 2 , 
giving |u 2 | = |w 3 | = 1. Thus, let us consider 

Wj ( m ) = — sin 7+ (—1 ) J ie ltPrn cos 7 

with j = 1, 2 and concentrate on the final equation sin a — Wj(m ) cos a = e lm Wj(m). 

Clearly, if | cosa| = 1, then the matrix is simple enough to have a biunimodular 

Hence, we assume |cosa| < 1 and claim that ^ has a desired 
as long as we can find toq £ R and jo € { 1 , 2 } such that 


vector, see Remark 6.1 
solution x, y 1 to. £ 


Wj o {m 0 )\ = 


Slllft 


gtrno 


COS OL 


(9) 


Slllft 


Indeed, after establishing a solution to ©, we can find x$ £ R such that 

etx °Wj 0 (mo) = 

3 ' e im ° + cos a 

so sina — e lx °Wj 0 (mo) cosa = e m e lx °Wj 0 (mo), leading to the conclusion that mo, jo and xo 
solve © by setting (e tx ,e ly ) = e“°(l, (—l)/°ie lVm o). 

Thus, we are left with solving ©• Here, it is important to notice that the biunimodular 
vectors (1 ,ie l(Pm ) and (1,— ie l(p7rl ) are orthogonal. Thus, if we define fj(m) = |u 7 (to )| 2 for 
j = 1,2 and to £ [0, 27r], then fi(m) +/ 2 (m) = 2||(— sin 7 , cos 7 )||§ = 2. Moreover, there 
is an mi £ [0,27t] such that e lmi = z , and then ip mi = 0 by ©, so f\{m\) = / 2 (toi) = 1. 

Finally, let us notice that for g(m) = 

SU P 9= 1 —| cosa| J > L 


' l +cos a 


one has inf 5 = ( ^TosLl ) ^ 1 and 
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This allows to finish the proof. Indeed, if fi > g and / 2 > g, then 1 > g contradicting 
sup g > 1. Similarly, if f\ < g and /2 < g, then 1 < g contradicting inf g < 1. Thus, 
there is an m 2 £ [0, 27t] such that < 7 ( 7712 ) is between f\{m. 2 ) and /2 (t 712)- However, g can not 
be strictly between / 4 and /2 on [0, 27 t], because fi(m\) = fiimi)- Therefore, there is an 
mo £ [0, 27 t] such that either < 7 ( 771 . 0 ) = /i(i77o) or g(mo) = / 2 (mo). □ 

Remark 6.4. Since for any vector R £ C n with || 1?.||2 = 1 the average / T „ \(R,u)\ 2 du is 
equal to one, it is easy to see, that for such R there is a u £ T n , such that |(1?, u)\ = 1. 

For two orthonormal vectors R\,R 2 £ C” finding a vector u £ T" such that |(2?i,m)| = 


\(R 2 ,u)\ = 1 seems to be a much harder task, and Theorem 6.3 is equivalent to accomplishing 


this task with n = 3. Currently, the only arguments extending this property for n > 4 rely 


on Theorem 3.f For orthonormal R\,R 2 £ £ 1 (Z) C ( 2 (fL) the problem is open. 


As in the 17(2) case, the existence of biunimodular vectors for all members of U (3) allows 
to write a formula for an arbitrary matrix in U( 3). 

Corollary 6.5. A £ 17(3) if and only if 

0 

a(l + z ) 

6(1 - z) 


A =\ 


0 

0 

>< 


'1 1 1' 


0 A 2 0 


1 UJ (jJ 


1 

CO 

0 

0 

_1 


1 U ! UJ 



0 


ac(l — z) 
6c(l + z) 



A 

1 

1 


A 

0 

o' 


1 

UJ 

uJ 


0 

a 4 

0 


1 

LU 

UJ 


0 

0 

A5_ 


for some parameters a, 6, c, z, Xj £ T, j = 1,..., 5 and ui = e 7 . 

Proof. Let A £ U( 3). By Theorems |6.3| and 4.1 we have that A £ V • Fix(l) • 2?! 


Moreover, in this case Fix(l) = F 3 ■ Fix((l, 0,0)) • F£, where F :i is the Fourier transform and 


Fix((l,0,0))= < 

h 

'20 o' 

0 a(l + z ) ac(l — z) 

: a, 6, c, z £ T 


l 2 

0 6(1 — z) 6c(l + z) 

J 


due to the description of 17(2) provided in Corollary 




A 

1 

1 " 

5.2 

Since F :i = X 

1 

UJ 

UJ 


1 

UJ 

UJ 


□ 


we can finish the proof. 

Moreover, as a byproduct from the proof of Theorem |6.3| we obtain the following descrip¬ 
tion of 17(3). 

Corollary 6 . 6 . 


17(3) = {diag( y \i,\2,\3,)XpYf t X 1 diag(l, A 4 , A 5 ) : a,/3,7 £ [0, 27t]; z, Xj £ T, j = 1,..., 5} 
= {diag(\i, \ 2 , \ 3 )Z l 3 Y*X 1 diag(l, A 4 , A 5 ) : a,/3,7 £ [0, 27 t]; z, Xj £ T, j = 1,..., 5}, 


where for a given ip £ [ 0 , 2 it] 



A 

0 

0 


COS pj 

0 

— sin p 


cos p 

— sin p 

o' 

x v = 

0 

cos tp 

— sirup 

11 

N Cl 

0 

z 

0 

, Ztp — 

sin p 

Q,OSp 

0 


0 

sirup 

cos <p 


sin p 

0 

COS(/9 


0 

0 

1 


Proof. This follows from Observation |A.3| used in the proof of Theorem |6.3| and the 
matrix factorization: 


'1 

0 

0 


cos a 

0 

— sin a 


a 

0 

0 

0 

cos 0 

— sin/3 


0 

z 

0 


0 

cos 7 

sin 7 

0 

sin/3 

cos /3 


sin a 

0 

cos a 


0 

sin 7 

— cos 7 


that leads to the provided description of U (3) via easy matrix manipulations. □ 

The above description of 17(3) can be compared with the standard Euler-Tait-Bryan 
decomposition of 0(3): 

0(3) = {X P Y*X 1 : a, 0 ,7 £ [0, 2t t]-,z = ±1} = {ZpY^ : a,/?, 7 £ [0, 2t r];*r = ±1}. 


































13 


7 U( 4) and beyond 


The study conducted in the two previous sections allows to surmise that, excluding some 
trivial cases, the number of biunimodular vectors for A £ U(n) is finite when n < 3. The 
classic example of the Fourier transform F 4 indicates a crucial change in the amount of 
biunimodular vectors, when n > 4. To understand this phenomenon it is enough to notice 
that (1, 0,1,0) and (0,1,0, —1) are eigenvectors of F 4 . Thus, (1, z, 1, — z) is biunimodular for 
F 4 for all z £ T. An investigation of the remaining biunimodular vectors for F 4 , that are 
given as (1, z, —1, z), allows to see a pattern that guarantees a continuum of biunimodular 
vectors. 

In order to view this pattern, for u £ C n we define its support as supp(u) = {k £ 
{ 1 , 2 ,..., n} : Uk 7 ^ 0 } and consider the following. 

Observation 7.1. Let A £ U{n). If there are u,w £ C" such that each of the vectors 
|u|, |u>|, \Au\, |Aui| is equal to 1 on its support and both disjoint unions supp{u) U supp{w) = 
supp(Au) U supp(Aui) cover {1,2,..., n}, then u + zw is biunimodular for A for all z £ T. 


The above observation allows to construct matrices in U(n) with all nonzero entries 
having a continuum of biunimodular vectors, whenever n > 4. It also allows for a further 
appreciation of Haagerup’s result stated in Theorem |1.2| It is not clear, however, whether 
for all unitary matrices with a continuum of biunimodular vectors one can find vectors u and 
w as in Observation 7.1 Even if this would be true, pointing to a conclusion that for any 


unitary matrix the set of its biunimodular vectors is a (finite) union of tori, it still would not 
tell, whether the set can be empty or not, underscoring the importance of Theorem |3.4| 

In the next section, we provide an alternative way of building U(2 n ), that is based solely 
on the biunimodular structure of t/(2). After that, we conduct the study of biunimodular 
vectors in the general U{n) case from the Lie groups perspective. Finally, we present a 
numerical treatment of the biunimodular problem with certain applications to the classical 
Fourier case. 


8 U{ 2 n ) 


In this section we show, that a further generalization of the biunimodular problem allows to 
build all matrices in U(2n) from matrices in U(n). 

In order to achieve this goal we employ the formula for U{ 2) given in Section [ 5 } 


fl 

a 

o' 


\2 

0 

b 



1 + z 

1 - z 


1 - z 
l + z 



1 O' 


0 c 


: a, b, c, z £ T 


( 10 ) 


and replace the unimodular numbers a,b,c,z£ T by unitary matrices A,B,C,Z £ U ( n ). 
Theorem 8.1. For every n £ N 


fl 

A 

O' 

'I+ z 

I - Z 

7 

o' 

12 

0 

B 

I - Z 

I+ z 

0 

c 


: A, B,C,Z £ U(n) 


( 11 ) 


where I stands for the identity matrix in U (n). 

The above theorem yields immediately a map from T 4 onto U(2 n ) via an obvious in¬ 
ductive procedure starting from the trivial map from T onto 17(1). In the first step of this 
procedure we recover (10). In the next step we get a formula for 17(4) that is enclosed in 
Appendix A. 

Before we prove the theorem, we need a slim primary on block matrices. Let M(n ) denote 
the set of all n x n matrices with complex entries. An arbitrary matrix in M(2n) can be 

^ ^ with A, B, C, D £ M(n). 

O L) 


written in a block form as 
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'A B 


A' B' 


AA' + BC' AB' + BD' 


A B 

C D 


C' D' 

— 

CA' + DC' CB' + DD' 

and 

C D 


Clearly, we have that 

. Moreover, every 2 nx n matrix with complex entries can be written as 
A', Y £ M(n) and 


C* 

D* 


with 


A 

B 

X 


AX + BY' 

C 

D 

Y 


CX + DY 


Lemma 8.2. Let -4= 

V 2 

U(n). 

Proof. Let t/ = -L 


A B 
C D 


£ U (2 n) with A £ U(n) and B,C, D £ M(n). Then B,C, D £ 


A B 
C D 


and I £ U{n) be the identity matrix. Since UU* = 


I 0 
0 / 


get that BB* = I. Considering U*U = 


I 0 
0 I 


yields that C*C = I and D*D = I. 


we 


□ 


Proof of Theorem \8.1\ We can recognize the task as a special sort of a biunimodular problem. 
For an arbitrary unitary matrix U £ U(2n ) we need to find a unitary matrix C £ U(n) such 
I 


that U 


C* 


is “unimodular”, in the sense that 


U 


with A,B,C £ U ( n ). 


( 12 ) 


(The general biunimodular problem of finding A,B,C,D £ U ( n ) such that U 


D* 


A 

C* 


B 


collapses to the latter since U 


D* 

C* 


I 

C*D 


AD 

BD 


.) Therefore, the goal is to 


iff U 

repeat the arguments used in Section [5 carefully replacing the unimodular numbers in T by 
unitary matrices in U(n). 

- - - - -j 

I 

I _ 7 

I ~ I 


'A* 

0 ' 

U 

7 

0 ' 

0 

B* 

0 

c* 


We start by noting that solving (|12|) amounts to saying that 

= \v £ U{2n) : V 


r/i 


\A* 0 


I 0 


I 

i 

, i.e. that V = 

0 B* 

U 

0 c* 

belongs to Fix 

I 


The description of the group Fix 


is again related to Fix 


I 



7 

o' 

0 



0 

A 


since F = ^ 


I I 
I -I 


£ U(2n) and F 


7' 

. 1 

7' 

0 

V2 

I 


being in turn equivalent to V = F 
In this way we get that 




I 0 
0 Z 


F = 1 
2 


Thus, V £ Fix 


: Z £ U (n) 

iff FFF £ Fix 


I + Z 
I - Z 


I - z 

I + Z 


with some Z £ U ( n ). 


A 

o' 


0 

B 



I + Z I-Z 
I-Z I + Z 



7 

o' 


0 

c 


for some A,B,C £ U(n) iff (121 holds. 


In order to solve the biunimodular problem (12) we write U = 


X Y 

X I Y' 


with A, y, A', Y' £ 


M(n ) and apply the singular value decomposition (SVD) to X and Y. 

Recall that for an arbitrary matrix M £ M ( n) there are unitary matrices W, W' £ U (n) 
and a diagonal matrix with non-negative entries E £ M{n ) such that M = W'YW*. This 
leads to the polar decomposition M = SmUm with Sm = W'Y>(W')* - a positive semi- 
definite matrix and a unitary matrix Um — W'W* £ U(n). 
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In short, we use polar decomposition A' = S X U X , Y = SyUy to write 

U = 

and claim that 

U 


SxUx SyUy 
X 1 Y' 


I 


A 

}U*U X 


B 


that A={S X + iS Y )U x and AA* = S 2 X + Sf + i(S Y S x - S X S Y ). Since UU* = 


with A,B£ U(n) as desired, so C* = iU Y Ux solves problem (12). To see this, we check 

7 O' 

0 I we 

get that S x + S Y = /, so S x and S Y commute. For positive semi-definite matrices Sx and 
Sy this means that they commute as well and, therefore, A £ U(n). 

To finish the proof we make the final observation that A £ U(n) already implies that 


B £ U{n). Indeed, we can extend the matrix ^ 

U[2n) and see that =U 
assures that B £ U(n). 


I C 
C* -I 


' I 

C' 

1 

A 

X"' 

C* 

- 1 

\/2 

B 

Y " 


to a unitary matrix ^ 

£ U(2n). Since A is unitary, Lemma 


8.2 


Therefore, we have proved that every U £ U(2n) can be written in the way prescribed by 

I 


(11). On the other hand, we have already checked that F 
U (2 n) for an arbitrary Z £ U (n), so formula © holds. 


0 


F=| 


I+Z I-Z 
I-Z I + Z 


£ 

□ 


9 U{n) 

We are interested in finding the biunimodular vectors associated to a matrix A £ U(n) for 


any n £ N. As it was pointed out in Remark 3.2 we can concentrate on biunimodular vectors 
with the leading entry equal to 1. Let us denote the collection of such vectors as Bi(A), so 

Bi(A) = {(v 1 ; v 2 ,.. .,Vn) £T n n A* (TP) : Wl = 1}. 

Moreover, let us make the following formal definition: 

B n = {A £ U(n) : Bi(A) ^ 0}. 

Clearly, B n is the set of n x n unitary matrices that possess biunimodular vectors and 


Theorem 3.4 asserts that B n = U(n). While it would be hard to reproduce this result, one 
can check what information about B n and Bi(A) can be drawn by using basic tools of Lie 
groups theory. 

9.1 The main correspondence 

The main tool for understanding B n and Bi(A) will be the factorization conducted in 
Proposition 3.1 As we noticed in Theorem |4.1[ a) the proposition implies that A £ B n 
iff A £ V ■ Fix(l) • T>i, where V C U{n) denotes the subgroup of diagonal matrices with 
diagonal entries in T, T>i C T> is the subgroup of matrices with the first diagonal entry equal 
to 1 and Fix(l) = {S £ U(n) : S 1 = 1} with 1 = (1,1,..., 1) S C n is also a subgroup of 
U{n). This allows to consider the following map 

$ : V x Fix(l) xPi-> U(n ) 

${D 1 ,S,D 2 ) = D 1 SD 2 
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containing the entire information concerning the existence of biunimodular vectors. 

Indeed, the matrices possessing biunimodular vectors are given as the image of $, what 
already provides some information on B n . 

Proposition 9.1. B„ = V • Fix(l) ■ V i = x Fix(l) x V i). In particular, B n is a closed, 
pathwise connected subset ofU(n). 

Proof. As we mentioned before, the equality B n = T> • Fix(l) ■ V\ follows immediately 
from Theorem |4.1[ a). Thus, by definition, B n = $(2? x Fix(l) x T>i). Since the subgroups 
V, Fix(l), V i are compact and pathwise connected, B n , which is the continuous image of 
their cartesian product, has the same properties. □ 

Moreover, the preimage <h -1 (A) bijectively corresponds to Bi(A). 

Proposition 9.2. There is a continuous bijection between Bi{A) and <b _1 (A). 


Proof. With the notation that D u is a diagonal matrix with entries on the diagonal 
given by u £ C n , the bijection can be described as follows. The preimage $ _1 (A) consists of 
all triplets (Hi, S', D 2 ) £Vx Fix(l) x T>i such that A = D\SD 2 . In light of Proposition 3.1 
this means that D 2 = Dy with v £ Bi(A), D x = Da v , and S = Df^ADif 1 = D^AD V . 
Therefore, the projection (D\, S, D 2 ) -A D 2 and the map D u -A u yield the continuous 
bijection from $ _1 (A) onto Bi(A), that is given by (Hi, S, D u ) -Au. □ 


Remark 9.3. U(n) has (real) dimension n 2 , whereas Fix(l) = U(n — 1) is of dimension 
(n — l) 2 . V has dimension n and T> i has dimension n — 1, hence $ is a mapping between 
two manifolds of (real) dimension n 2 . Clearly, $ is smooth. 


9.2 The structure of B n 

In this subsection we show, that for every n £ N the set B n contains a nonempty open subset 
of U(n). This is done by calculating the Jacobian of the map $ : Hx Fix(l) x V\ —>■ U (n) and 
proving that at a certain point it has a full rank. This allows to apply the inverse function 
theorem to draw the above conclusion. 

In order to proceed with the calculations we provide some basic facts from the theory of 
Lie groups and establish the necessary notation. 

We consider U(n) and its subgroups as (real) submanifolds of the space gl(n, C) of ar¬ 
bitrary matrices. Given a closed subgroup H C GL(n, C) and g € H, the tangent space 
T g H is the space of all tangent vectors of curves in H going through g , conveniently viewed 
as elements of the ambient space gI(n,C) containing H. The tangent space T e H is a Lie- 
subalgebra of g[(?r, C), i.e. it is closed under the Lie bracket [X, Y] = XY — YX, and it is 
referred to as the Lie algebra () of H. Here, e is the neutral element, that is, the identity 
matrix that we shall denote by I. An alternative description of 1) is obtained via the matrix 
exponential, defined by 

k—0 

It is easy to verify that (f exp(tX)) |j =0 = X, and fi consists of all matrices A' such 
that exp(tX) £ H , for all t £ R. 

Since multiplication (on the left or on the right) with a group element g is a diffeomor- 
phisrn on H, that is a matrix group, one obtains that 

TgH = g ■ b = {gX : X £ f)} = f) • g = {Xg : X £ b}. 

We let u,0,0i,f denote the Lie algebras of U{n), V, V\, Fix(l), respectively. Then 

u= {X £ 0 [(n,C) : A* = -A} , 
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5 consists of all diagonal matrices with purely imaginary entries, and Iq C 0 is the subspace 
of matrices with vanishing upper left corner. Differentiating the equality exp(tJT)l = 1 at 
t = 0 yields the following description of f: 

f = {Ieu:Il = 0 }, 

i.e. the elements of f are characterized by the fact that the sum over rows are zero. We also 
note the simple fact that 5 fl f = { 0 }. 

Finally, recall the definition of Jacobians. If F : X —> Y is a smooth map between 
manifolds X and Y. the Jacobian of F at x £ X is the linear map dF(x) : T x X —> T F ( X \Y 
defined as follows: For any v € T x X pick a curve 7 : (—e, e) — > X with 7 ( 0 ) = x, 7 ( 0 ) = v. 
Define the curve 7 = F 07 : (—e, e) —> Y (by definition 7 ( 0 ) = F(x)). Then dF(x)(v) = 7 ( 0 ). 
To calculate the Jacobian of $ : V x Fix(l) x T>i —► U(n) we shall use the following 

Lemma 9.4. Let Hi, H 2 be closed subgroups ofU(n) and m : H 1 x H 2 —> U(n), m(g,h) = 
g ■ h. Then, for all g,h £ U(n), the Jacobian of m at ( g , h) is given by 

dm(g,h)(X,Y) = gY + Xh , VX eT g H u Y eT h H 2 . 

Proof. Let 7 , : (—e, e) —> Hi denote smooth curves with 

7i(0) = <7)7i(0) = X and 72 ( 0 ) = h, 72 ( 0 ) = Y . 

By a standard reasoning, the product rule for scalar-valued functions easily extends to 
matrix-valued functions implying that 

d -^mm {0) = 7 ^ 0 ) 72 ( 0 )+ 71 ( 0 ) 72 ( 0 ) = g Y+xh. 

□ 

Using this lemma and the chain rule, we obtain a rather transparent description of the 
Jacobian d$(Di, S, D 2 ): 

Lemma 9.5. Let {Di,S, D 2 ) GVx Fix(l) x V 1 . Then 

d${Di, S, D 2 ) : T Dl T> x Tg(Fix(l)) x T D 2 V 1 —y T Dl sD 2 U(n ) 


acts via 

d$(D u S, D 2 )(X, Y, Z) = Di(SZ + YD 2 ) + XSD 2 
Proof. Since we can factor 


$ : (Di,S,D 2 ) h> ( D 1 ,m(S,D 2 )) m{Di,m{S,D 2 )) , 

applying the previous lemma twice and using the chain rule yields 

d$(D 1 ,S, D 2 )(X, Y, Z) = [dm(D 1 ,SD 2 )o(I,dm(S,D 2 ))](X,Y,Z) 

= dm(D 1 ,SD 2 ){X,SZ + YD 2 ) = D^SZ + YD 2 ) + XSD 2 . 


□ 

Recall that we are particularly interested in the rank of that is, the real dimension of 
the image of dd>. The following lemma translates this to a more tractable problem in linear 
algebra: 

Lemma 9.6. Let S £ Fix(l) be given, and m £ N. Then the following are equivalent: 

(a) For all (any) ( Di,D 2 ) £ V x V 1 the Jacobian d<&(D\, S, D 2 ) has rank m. 

(b) S, I) has rank m. 
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(c) The linear map 


C s :dx f x 5i ->• u , (X,Y,Z)m-X + Y + SZS* 


has rank m. 

Proof. For the equivalence (a) (b) we observe that 4>(CT?i, S, D 2 C') = C$(Di,S, D 2 )C' 

for all C GV, and C' G V\. Thus, by the chain rule, d<&{CD\, S, D 2 C') = Cd${D\,S, D 2 )C ', 
what yields the equivalence. 

To see that (b) <t=> (c) we note, that by the previous lemma, d<E>(7, S, I)(X,Y' , Z) = 
XS + Y' + SZ, for all X G 0, Y' G Ts(Fix(l)), Z G hi. Since S is invertible, the rank of 
d$(I, S, I) equals the rank of Rs* o d$(7, S , I), where Rs* denotes right multiplication with 
S*. Clearly, Rs* od$(7, S, 7) (A, Y', Z) = X + Y'S* + SZS* and this map has the same rank 
as the map Cs from (c), since Y = Y'S* G f iff Y' G Tg(Fix(l)). □ 

We next translate the condition from part (c) of the previous lemma to a simpler one. 
In the following lemma and its proof, the real and imaginary parts of complex vectors and 
matrices are computed componentwise. In particular, given a complex-valued matrix U, 
Im(f7) is obtained by taking the imaginary parts of the individual entries of U and the rank 
of Im({7) is its rank as a mapping on R". Also note that throughout the proof of the lemma, 
“linear” means “R-linear”; and matrices are identified with the linear map they induce. 

Lemma 9.7. Let S G Fix(l), then 

rank(d4>(7, S, I)) = (n — l) 2 + n + rank(Im(S)) . (13) 


Proof. Let Cs : h x f x hi —> u denote the linear map from Lemma [976} Let n = h + f C u. 
Then h Cl f = {0} implies that the restriction of Cs to h x f is an isomorphism onto n. This 
allows for the first step towards (131 by considering the quotient map 


C s '■ hi u/n , Z ^ SZS* + i 


By Lemma |9.6| and a rank formula for the quotient map we obtain that 

rank(d$(7, S, 7)) = rank(C's) = rank(C l | 8X f) + rank(Cs) = (n — l) 2 + n + rank(Cs) . (14) 

In order to prove that rank(C's) = rank(Im(S 1 )), we need to free ourselves from the 
quotient space u/n. For this purpose, let us consider the map T : u —> C n , T(X) = XI. We 
claim that 

T(u) = V :={xGC n : Re(x, 1) = 0} . 

Indeed, the kernel of T is precisely f, which has real codimension 2n — 1 = diniR(P). Thus 
it suffices to prove that T(u) C V. This, however, follows from the fact that A' = —A*, for 
all X G u, since 


Re(T(A), 1) = Re(Al,l) =Re(l,X*l) =Re(—Al,l) =—Re(T(A),l) . 


We next consider the linear map J : u H R", J(X ) = Re(T(A)). The kernel of J is 
then given by n, as can be seen by yet another dimension argument. The observation that 
T(u) = V yields that 


J(u) = {Re(u) : v G T(u)} = {in G R n : (w, 1) = 0}, 


so J(u) is a subspace of R" of dimension n — 1. Since the real codimension of n is also 
n — 1, it is enough to show that n C ker(J). For the latter, note that J(X + Y) = J(X) = 
Re(xi,... ,x n ), for all Y G f and X = diag(xi,..., x n ) G h. Recall that 1 6 5 entails in 
addition Xk G iR, whence J(X + Y) = 0, establishing that ker(J) = n. 
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This allows to conclude that J induces an injective linear map J : u/n — > R". Hence, if 
we define Q = J ° Cs, then we obtain that rank(Cs) = rank(Q). 

Clearly, Q : hi —> R" and this map is much easier to treat than Cs- Indeed, given any 
Z £ 5i, we employ 5*1 = 1 to compute that 


Q(Z) = Re(SZS*l) = Re(SZl) = -Im(S)z , 


where z = (0, Z 2 , ■ ■ ■ , z n ) £ R" is related to Z via Z = diag(0, iz 2 , ■ ■ -, iz n ). This shows that 
rank(Q) = rank(Im(5)|w 0 ), where Wo C R ra is the subspace of vectors with first entry equal 
to zero. Moreover, Wo+R-1 = R” and Im(5)l = 0, since 51 = 1. Hence rank(Im(5)|w 0 ) = 
rank(Im(5)). 

Thus, we have established that rank(Cs) = rank(Im(5)), so equation (131 follows from 

©■ 

□ 

The above lemma leads to the fundamental fact, that for all n £ N the maximal rank of 
dd* equals the dimension of U(n). 

Lemma 9.8. For all n £ N there is an S £ Fix(l) C U{n) with rank(d<l>(/, 5, 1)) = n 2 . 


Proof. Let Wg C R n denote the orthogonal complement of R • 1 in R", and Win C C n 
denote the orthogonal complement of C ■ 1 in C n . In order to exhibit 5 £ Fix(l) with 
rank(d$(/, 5, /)) = n 2 let us pick S' £ Fix(l) n O(n), where 0(n ) stands for the orthogonal 
group. This entails 5 , (Wr) = Wr, so we can define a unitary mapping 5 on Wc ® C • 1 by 


S(v + al) = iS'v + al , a £ C, v £ Wc ■ 


Then 5 £ Fix(l), with Im(5)|w R = 5'|w R . Since 5'(Wr) = IFr, the restriction 5 'Iwk bas 
rank n — 1. Moreover, Wr © R • 1 = R” and Im(5)l = 0, so 


rank(Im(5)) = rank(Im(5)|u4) = rank(5 , |w K ) = n — 1. 


Thus, Lemma 9.7 implies that rank(d<l>(/, 5,/)) = n 2 


□ 


With the above lemma we can draw the mentioned conclusion regarding the structure of 


B„ 


Theorem 9.9. For every n £ N the set B n C U(n) has a nonempty open interior. 


Proof. By Propositi on |9 .1 B 
of dimension n 2 . Lemma 


9.8 


= d>(P x Fix(l) x T>i) and $ is a map between groups 
guarantees an existence of a point in V x Fix(l) x T >\, where 
the Jacobian of <F has the full rank. Therefore, by the inverse function theorem, $ is a 
diffeomorphism around this point. □ 


Corollary 9.10. B n = U(n) if and only if B n ■ B n C B n . 

Proof. Clearly, A £ B n iff A~ x £ B n . Therefore, B n is a symmetric subset of U(n) 
that, by the above theorem, has a nonempty open interior. If, in addition, it is closed under 
multiplication, then B n must be an open subgroup. Since U(n) is connected, it follows that 
B n = U(n). The converse is clear. □ 

The results of this subsection can be compared with a theorem of De Vos and De 
Baerdemacker based on Hurwitz parametrization and stating that ( B n ) n ~ 1 = U(n). 


9.3 The structure of Bi(A) 

Regarding the structure of Bi(H) it is crucial to recall Conjecture [T t] In light of the hypoth¬ 
esis imposed by Bjorck and Saffari, the main issue is to establish the cardinality of Bi(H). 
Therefore, there is a merit in proving the following. 
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Proposition 9.11. For almost all A £ U(n ) the set Bi(A) is finite. 

Proof. This follows from Proposition [972] and Sard’s theorem. By Proposition [972] we can 
replace Bi(A) by $ _1 (A). Clearly, we can concentrate on the case $ _1 (A) 0. Recall that 

the regular points of a smooth map are those for which the rank of the Jacobian is full; any 
point that is not regular is called critical. A regular value is a value, such that all the points 
in the preimage are regular. If A is a regular value of <F, then $ -1 (A) must be a discrete 
subset of the compact space V x Fix(l) x 'D 1 , hence finite. Thus, if $ -1 (A) is infinite, then 
A must be the image of a critical point of <b. However, Sard’s theorem states that the images 
of the critical points constitute a set of the Haar measure zero, what concludes the proof. 

□ 

In order to achieve a deeper understanding of the set Bi(A) for a given A £ U(n ), we 
consider the phasing manifold of A 


Ad A = 'DA'Di 


introduced in and the stabilizer group T>a = {(E\,E 2 ) £ V x V i : E\AE 2 = A}. 
By Proposition 9.2 the next result establishes a bijection 

Bi(A) —> ( Ma H Fix(l)) x T>a , 


that can be explained in the following way. Intuitively, finding all biunimodular vectors for 
A amounts to finding all S £ Fix(l) for which there exists a pair (D\,D 2 ) £ V x V\ such 
that A = D\SD 2 (i.e. S = Df 1 ADf 1 £ AAa) and, after that, finding other such pairs for 
S by considering the set {(D\Ei, D 2 E 2 ) : (£j . E 2 ) £ T>a}- We shall prove that this simple 
scheme yields all members of Bi(A) in a one-to one fashion. 

Theorem 9.12. There is a bijection W : (AAa H Fix(l)) x T>a —> <I > “ 1 (A). 

Proof. We can pick a map ip : AAa H Fix(l) —> V x V i, ip(S) = (ip\(S), ip 2 (S)), 
such that ipi(S) Sip 2 (S) = A. Indeed, by definition, for every S £ AAa D Fix(l) there is 
(. D il D 2 ) £ V x Vi, such that S = DiADo , thus we can pick one of such pairs and set 
(MS),MS)) = (DffiDf 1 ). 

Then, the map U/ : (AAa f~l Fix(l)) x T>a —t $ 1 (A) is given by 
T (5, (E 1: E 2 )) = (Eiipi(S), S, E 2 ip 2 (S)) ■ 

To see that the image of is contained in $ -1 (A), it is enough to note that 

$(T(5', (e 1 ,e 2 ))) = e 1 ms)sms)e 2 = e 1 ae 2 = a. 


The map is one-to-one, since T ( S , (Ei,E 2 )) = T (S' , (E(,E 2 )) implies that S = S' and this, 
in turn, implies that Ei = E[, for i = 1,2. To check that the map is onto, we take an 
arbitrary element (Gi,S, G 2 ) £ $ -1 (A) and easily verify that T (S,(Ei,E 2 )) = (Gi,#, G 2 ) 
with Ei = Gi (MS)) -1 , for i = 1,2. Since G 1 (MS))~ 1 A(MS))~ 1 G 2 = G X SG 2 = A, we 
get that (E\,E 2 ) £ T>a and end the proof. □ 


Remark 9.13. Since in the above theorem it is hard to guarantee the continuity of ip, we 
can not conclude that the bijection is continuous. We mention this, because otherwise, by 
Proposition 9.2 the set Bi(A) would be homeomorphic to (AAa HFix(l)) x T>a- 

With Theorem 9.12| we can concentrate on the set (AIa H Fix(l)) x T>a, in order to 
investigate the cardinality of Bi(A). For that, we need a better understanding of the phasing 
manifold AAa , for a given A £ U(n). This can be obtained by defining the mapping 


$A : P x T>i —> U(n) , a(Di, D 2 ) — DiAD 2 


and noting that AAa — ^a(P x T>i). This setup allows for making some useful observations. 
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Lemma 9.14. For every A £ U(n ) the set A 4a is the orbit of A under the left action of 
VxD i on U(n ) given by ((Hi, D 2 ), A) 1 —> D\AD 2 - In particular, Ma is a closed submanifold 
ofU(n) with Ma — (V x V\) /Da, so dim(AlA) + dim (Da) = 2n — 1. 

Proof. Note that the map ((D\, D 2 ), A) 1 —>■ D 1 AD 2 indeed defines an action of the 
product group, since the diagonal group is commutative; and clearly Ma is the orbit of A. 
The remaining statements follow from this: is just the quotient map with respect to this 

action, what implies that the differential of 4 >a has constant rank; and this in turn allows 
to define a manifold structure on Ma- In addition, Ma is compact, hence closed, and *I>a 
induces a continuous bijection (V x T>i) I'D a —> Ma between compact Hausdorff spaces, 
which then has to be a homeomorphism. Therefore, dim(AfA) +dim(I?A) = dim(2? x Di) = 
2 n — 1. □ 

Corollary 9.15. dim(AlA) = 2n — 1 dim(I?A) = 0 <!==> Da is finite. 

Proof. The first equivalence follows immediately from the dimension formula given in the 
above lemma. In the second equivalence one implication is obvious. For the other implication 
we note that, if dim(2?A) = 0 then Da must be a discrete subgroup of the compact group 
2? x 2?i, so it must be finite. □ 


At this stage we can provide the following explanation on the cardinality |Bi(A)| of Bi(A). 

Obviously, for n = 1 we have |Bi(A)| = 1. When n = 2, then |Bi(A)| equals to 2 or to the 
continuum c. For n = 3 we have only found examples where |Bi(A)| equals to 4, 5, 6 or c. By 
Theorem 9.12 |Bi(A)| is equal to the cardinality of the set (Ma H Fix(l)) x Da- Therefore, 
to see the general picture, we employ the above corollary together with Theorem |3.4| and get 
the following: 


1. If dim(AlA) < 2n — 1 then |Bi(A)| = c. 

2. If dim(AlA) = 2n — 1 then 

(a) |Bi(A)| e N or 

(b) Bi(A) is infinite but countable (no examples) or 

(c) Bi(A) is uncountable. 

Remark 9.16. Regarding the case clim(A4A) < 2n —1, it is clear that for the identity matrix 
I (and any diagonal matrix) we have that Mi = D, so dim(A4/) = n. Moreover, there are 
examples of matrices A £ U(n) with dim(AlA) in the range from n to 2n — 1. 

We would like to end this section by getting closer to understanding the classical Fourier 
case. As it turns out, for the Fourier matrix F n the stabilizer Dp n is trivial. 

Lemma 9.17. For every n £ N we have that Dp n = (I, I). 

Proof. Let (Hi, £ 2 ) £ Dp n , i.e. (Hi, £ 2 ) £ D x V\ and EiF n E 2 = F n . Then 
£2 = F*E(F n . This means that £2 commutes with all translations (time shifts) and all 
modulations (frequency shifts). A result going back to Schwinger asserts that translations 
and modulations generate a basis for the space of all matrices 0 t(n, C)Q Therefore, £2 com¬ 
mutes with all matrices in this space, so £ 2 = ±1. Since E 2 £ D 1 , we get that £ 2 = /, so 
£1 = I as well. □ 

The above lemma together with Theorem |9.12| and Lemma |9.14| entail the following 
information on the set of all biunimodular vectors in the classical Fourier case. 


Theorem 9.18. For every n £ N there is a bijection between Bi(F n ) and Mp n (~l Fix(l), 
where M F n is the phasing manifold of dimension 2n — 1. 

4 If T is a time shift by one, and M is a similar frequency shift, then T k M l : k,l = 1 ,... ,n} form an 
orthonormal basis of £)l(n, C) (see [47]). 
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The above theorem indicates the resilience of Conjecture 0 Another approach towards 
this conjecture has been designed by Gabidulin. Unfortunately, even in the prime case, where 
some details are provided, he does not explain how to exclude that Bi(.F n ) is infinite but 
countable. Therefore, the only accessible confirmation of the conjecture in the prime case is 

m- 


10 Numerical Results 

10.1 A simple algorithm and its properties 

In the proof of Theorem |4.1[ b) we have shown that H-AHoo^i < n for every A £ U (n). As we 
pointed out in Remark |4.2[ biunimodular vectors for A are precisely those, where the norm 
|| A||oo— >-i == n is attained. Thus, if v is a biunimodular vector of A , then we can think of it 
as of a point on T", where the norm ||Av||i is maximized. 

In order to maximize the £} -norm it is convenient to notice, that for every v £ C n with 
\v\ < 1 (entrywise) and A £ U(n) we have that 

||Au||i = (Av, sign(Au)) = (v, A*(sign(Au))) < || A* (sign(Au))||i, (15) 


where “sign” is the standard signum function defined on the support of v £ C" by 

sign(v) = 

M 


so that v = sign(u)|u|. Since (15) holds for all v with |u| < 1, we can safely define sign(u) as 
zero on the complement of the support of v. 


Applying (151 twice yields 


II^Hli < ||A*(sign(Au))||i < ||A(sign(A*(sign(Au))))||i := ||A/||i, 


(16) 


with v' = sign(A*(sign(Ai>))). This leads to the following algorithm that with each step 
increases the ^ 1 -norm of A on consecutive vectors, and therefore, allows to search for the 
biunimodular vectors of A. 


1. Initialization: Take any vector Vq £T n . 

2. Step: Vj +1 = sign(A*(sign(AU i ))) for j = 0,1,.... 


Clearly, (16) assures that 


\\AVj\h < ||AU J+1 || 1 for i = 0,1 ,... 


(17) 


giving a way to stop the algorithm after reaching J such that ||AVj||i = n. Of course, for 
numerical algorithms we can only expect this to hold within a certain, predefined precision. 
Since we know that ||AVjj|i < n for all j £ N, we can set this precision to S > 0 and stop 
the algorithm after reaching a J such that 


||AUj||i > n- 6. 


(18) 


In this way we rediscover a phase retrieval method that has a pretty long history. To 
make it short, following Jaming 133 we recall that as early as 1932, Pauli was asking whether 
the information on the moduli of a wave function and its Fourier transform allows to recover 
the phase of the wave function (i.e. Having 1^1 and | can we recover ip?). In practice, the 
phase retrieval problem is a serious obstacle for analysing complex valued signals, since only 
partial information is given on |^| (see p2J [IS] by Candes et al. or [1] and for recent 
developments on this notorious problem). Nevertheless, under original Pauli constrains, a 
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method for recovering the phase of a signal has been provided in 1972 by Gerchberg and 
Saxton (see [32). Gerchberg-Saxton algorithm has several variants, and the one that is 
relevant to us takes the following form. Let \ijj\ = g and \Pip\ = h. In order to recover %p 
define 

Pghf = 3 sign 1 (J‘*(/isign 1 (J : '/)))), 


where 


signi(/) 


| 7 | on supp(/) 
1 off supp(/) 


and starting with an fo such that |/o| = g, run fj = P gh (fo) until it converges to a solu¬ 
tion This algorithm can be viewed as an alternating projection method for finding an 
intersection of two subsets of a Hilbert space, that was discovered by von Neumann in 1933 
(E2)j hr the case when the subsets are closed subspaces. This view is justified by writing 
P g h as Pq o Ph, where Pcf = gsign 1 (/) and Pnf = T*{h .sign 1 (J 7 /))) are projections on 
certain subsets of a Hilbert space. 

The method of alternating projections is a topic by itself (see [26]), that is well covered 
in the case of convex sets ([5]) and compact sets with a non-empty intersection ([21]). In our 
setting we are looking for an intersection of two compact non-convex sets, namely T n and 
A*(T n ), with a guarantee that T n n A*(T n ) ^ 0 provided by Theorem 

Clearly, signj : C n —> T n is a projection on the torus (in the sense that ||i> — sign 1 (u )||2 = 
dist(u, T n )), and A* o sign 1 oAisa projection on H*(T”). Thus, a successive application of 
P' A { V ) = sign 1 (H*(sign 1 (Hu))) yields the standard alternating projection method. 

Although we use 


3.4 


Pa{v) = sign(A*(sign(Au))), v £ C", A € U(n) 


instead of P' A , we shall quickly confine our search to the range where both these operators 
agree. Thus, we shall refer to this method as “alternating projection algorithm” anyway. 

After providing this backdrop, we can take on understanding the workings of the algo¬ 
rithm. Recall that we start with an arbitrary vector Vq £ T n and successively apply Pa, by 
setting Vj = P a Vq, until we reach a J such that (181 holds with a certain precision S > 0. 

The price for the simplicity of this algorithm is its randomness, that has to be properly 
interpreted]^] 

Since v £ T n is biunimodular for A iff ||Av||i = n, we would be satisfied with finding 
Vj such that S in (18) is close to zero. However, there is no guarantee that for a given 
starting vector Vq we will find J £ N so that S is small enough. This is already apparent 

'1 1 " 


in the case A = Fo = 


l 

V2 


1 -1 


with V 0 = (1,1). Indeed, here Vj = V 0 for all j £ N, 


so ||AH, ||i = y/2 <C 2. Therefore, in general, the algorithm has to be applied to multiple 
random starting vectors in T". 

Another aspect of the randomness is due to the discontinuity of the signum function. 
Clearly, neither P A nor P' A is continuous. In theory, the vector 1 = (1,1,..., 1) £ C n is a 
fix point of both P A and P A , when A = F n . In practice, both algorithms may converge to 
more stable solutions when applied to the starting vector 1 in the Fourier case. Of course, 
similar problems may occur in the general case as well. In order to avoid the discontinuity, 
we observe that for every A £ U(n ) and v £ T" such that 11 Az; 11! > n — S we have the 
following 

II 1 - I Av 111 2 = n - 2||Au||i + ||A V ||! = 2 (n - ||Au||i) < 28, 


5 Unfortunately, ip 1 is not unique (it depends on the choice of fo), so it is only a candidate for tp (see Subsec¬ 


tion 


10.31. 

6 Essentially, one has to run the algorithm with multiple starting vectors in T n to get satisfactory results (see 
Subsection 


10 . 21 . 
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so \Av\ > \ for 8 < |. Moreover, for such <5 we get 

||A*(sign(Au)) - v\\j = ||sign(Au) - Av\\\ = ||1 - |^4v|||i < 28, 

so |A*(sign(Ai>))| > | as well. This allows to draw the following conclusions. 

If the stopping condition (181 is satisfied, then || 1 — | AVj\ || 2 < V28- If the initial condition 

PVb||i>n-| (19) 


is satisfied, then for all j G N U {0} we have that 

. Vi := P J A V 0 = ( p' A yv 0 , 

• Vj G T n , \AVj\ > \ and |A*(sign(AV)))| > 5 , 


due to an easy induction argument. Moreover, assuming (19), it can be proved that 


\AV J+1 \\ 1 = \\AV j \\ 1 iSV j+1 = V j 


31 


9 H m j —>00 ||V)+i Ij'11oc 9. 

Where the former follows from (|T5| and (16), while the proof of the latter is a bit lenghtly, 


so we include it in Appendix B containing a detailed study of the algorithm. 

Since the sequence (||AV^||i)jgN is monotonic and bounded, it converges to a limit L < n. 
In Appendix B we show that under the initial condition (19) the cluster points of 


belong to the set Sa of fix points of P A . This allows to assure the convergence of Vj, in the 
case when Sa is finite. Since biunimodular vectors of A belong to Sa, we know that for some 
A the set Sa is uncountable. Thus, it it possible, that in some cases Vj does not converge. 

This brings us to the main point of the analysis of t he algorithm. If the stopping condition 
(18) holds with 8 < 4, then the initial condition (19) holds with Vo replaced by Vj. Thus, 


the stopping condition (18), with 6 numerically near to zero, gives Vj that is close to a fix 
point of Pa- If L := lim^oo ||AV)||i = n, then Vj is close to a biunimodular vector of A, 
however, it is not clear how close. If L < n, then Vj is close (in a similar fashion) to a fix 
point of Pa, that may be far away from biunimodular vectors. Moreover, it is impossible to 
establish numerically that L = n. Hence, we need to settle for 8-near biunimodular vectors 
given as 

Bi 5 (A) = {v G T" : ||Au||i > n ~ 5} 


and stress again, that some of these vectors can be far away from biunimodular vectors, even 
if 8 is nearly zero. 

In short, if the algorithm returns vectors that are “biunimodular” in practice, there is no 
guarantee that they are close to the theoretical biunimodular vectors. 

Despite this uncertainty, it turns out, that finding near biunimodular vectors for every 
unitary matrix is sufficient to conclude, that all such matrices have biunimodular vectors. 
In the following proposition we prove even more. 


Proposition 10.1. Let a < 1. If for every n G N and for every A G U(n) there is v G T n 
such that 

\\Av\\i > n - n a , (20) 

then every unitary matrix has a biunimodular vector. 

Proof. This follows from considering diagonal block matrices. If B G U(m) does not 
have biunimodular vectors, then = m — e with some e > 0. Therefore, for every 

N G N, the block matrix An with N copies of B on the diagonal 


An := 


B 


B 


B 
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belonging to U{n ) with n 
for Am, then we get that 


Nm, has the norm ||-Ajv||oo-^i = N(m 


e). Thus, if (20) holds 


n - n a < IIAatuHi <n- n-^, 


leading to a contradiction < n“ _1 , since n“ _1 —> 0 as N —> oo. □ 

In the next subsection we check the performance of the alternating projection algorithm 
for unitary matrices up to dimension 100. The final subsection contains results of the appli¬ 
cation of this algorithm to the Fourier case. 


10.2 Effectiveness of the algorithm 

In order to measure the performance of the alternating projection algorithm we recall, that 
for S > 0 and A £ U{n ) we have defined the set 

Bia(A) = {» e I" : ||Au||i > n - 5} 

of 5-near biunimodular vectors of A. For a given unitary matrix A and a prescribed parameter 
S the algorithm is supposed to return a member of Bi, 5 (A). The effectiveness of the algorithm 
can be measured by considering the set B n> s, that is, the set of these matrices in U(n), for 
which the algorithm returns the desired outcome. 

The aim of this subsection is to provide numerical evidence that the Haar measure of 
B nt s (that we denote by ix{B n , t 5 )) is very close to one. For this purpose we proceeded as 
follows: For each dimension under consideration, we drew a fixed number K of random 
unitary matrices with Haar measure used as underlying probability measure. To achieve 
this, we implemented the method described in m- We then ran the alternating projection 
algorithm on the matrix with randomly chosen starting vectors, in order to produce near 
biunimodular vectors. The search was terminated either if one such vector was found, or until 
a maximal number of starting points was exceeded. The parameters used in this procedure 
were as follows: 

• 6= 1CT 10 . 

• Dimensions: n = 3,5,10,25,50, 75,100. 

• Number of matrices: K = 10 4 ; 

• Maximal number of starting points: M = 10 3 ; 

• Maximal number of iterations: L = 10 4 . 

We implemented the described algorithms in MATLAB, and ran them on a stand-alone 
personal compute^ 

The results of our numerical experiment are documented in Table [1] Two facts seem 
particularly striking to us: First of all, the algorithm found a 5-near biunimodular vector 
for each of the 10 4 random matrices in each dimension considered. This translates to an 
estimate of fi(B nt s), which holds at least with very high likelihood: Our numerical experiment 
amounts to performing a Bernoulli experiment based on repeated independent samples of 
a Haar-distributed random matrix A in U{n ), and a positive outcome of the experiment 
(a near biunimodular vector being found) implies A £ B n Thus the probability of a 
positive outcome for a single event is given by p < p,( B n g )■ Assuming that fi(B n g) < 

1 — 10 -3 = 0.999, the probability of obtaining a positive result 10 4 times in a row will thus 
be < 0.999 10 « 4.5173 • 10 -5 . In other words, our experiments provide very convincing 
evidence that n(B nt g) > 0.999. 

Another striking phenomenon is illustrated by the 3rd and 4th column of Table |T] high¬ 
lighting the effectiveness of the search algorithm. For small to medium dimensions and in the 

' The CPU of the machine was a 6 core AMD FX 6100 processor with 8 MB level 1 cache, running at 3.3. 

GHz, together with 8 GB RAM. The operating system was SUSE linux 12.2 for 64-bit PCs. 
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Dimension 

Matrices in B n ^ 

starting points 
(average) 

starting points 
(maximal) 

total runtime 
(rounded) 

3 

10 4 

1.0663 

5 

1 min 38 sec 

5 

10 4 

1.1719 

7 

4 min 20 sec 

10 

To 4 

1.3472 

7 

9 min 59 sec 

25 

W 

1.6521 

10 

28 min 51 sec 

50 

W 

2.4106 

17 

1 h 34 min 

75 

W 

4.3361 

34 

5 h 6 min 

100 

10 4 

9.4730 

89 

16 h 40 min 


Table 1: Results of the numerical experiment 


average case, the algorithm needs only very few starting points to find a near biunimodular 
vector, but even the maximal number of such points stays well away from the maximal num¬ 
ber of 10 4 that was allowed. As the dimension increases, several effects conspire to increase 
the algorithm complexity: The rate of convergence in the alternating projection algorithm 
slows down, as witnessed by the increasing percentage of cases where the maximal number of 
iterations is reached before the (n — £)-threshold is passed (not shown in the table). Also, one 
may expect that the chances of randomly picking starting points for which the alternating 
projection algorithm converges sufficiently quickly also diminish with increasing dimension. 
Both effects increase the number of starting points needed in the search. In addition, the 
costs of applying the matrix-by-vector multiplication also can be expected to contribute to 
a nonlinear increase of the computational load. It seems that the overall result is an expo¬ 
nential increase of runtime: A log-linear fit for the runtimes for n > 10 revealed a behaviour 
like O(e 0 05n ). 

10.3 Biunimodular vectors for Fourier matrices 

The search for the biunimodular vectors for F n via the alternating projection algorithm is 
the simplest realization of the phase retrieval process under trivial Pauli constraints \ip\ = 
\T4>\ = 1 in the discrete Fourier case F = F n , if) G C n . 

When looking for elements of Bi(F„) for a Fourier matrix F n , it is advantageous to take 
into account various operations that preserve this set: 

• Circular shifts by k = 0,1,..., n — 1. I.e., if (uq, ..., u n - 1 ) G Bi(F„), then the vector 
(ufc,itfc_|_i,..., u n —\, Uq ,..., itfc- 1 ) is biunimodular. 

• Modulation by k = 0,1,..., n — 1: If (uq, ..., u n - 1 ) € Bi(Fji), then the vector¬ 
ed, ui exp(27rifc/n),..., u n -\ exp(27r ik{n — 1 )/n)) is biunimodular. 

• Dilation by k = 0,1,...,n— 1 coprime with n: If (uo, ■ ■ ■ ,u n -i) G Bi(F n ), then 
(yo, ■ ■ ■, Un- 1 ) with Dj = Ujk mod n is biunimodular. 

• Conjugation: If (u o,... ,it„_i) G Bi (F n ), then its conjugate (tZo,... ,u n _i) G Bi)! 7 ^). 

• Fourier transform: If u G Bi(F' n ), then F n u is biunimodular as well. 

These operations induce the action of a finite group G n of size 4 n 2 ip(n) on Bi(.F n ), where 
ip{n) is Euler’s totient function. 

We employed the alternating projection method to produce a list of near biunimodular 
vectors v G BR(.F n ), for 5 = ICE 7 . More precisely, we computed a list of orbit representatives 
of such vectors, by repeatedly running the alternating projection algorithm with random 
starting points. For all dimensions, we used 10 4 starting vectors, and at most 30 000 iterations 
before starting with a new random vector. Whenever the algorithm terminated in a d-near 
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biunimodular vector, the result was appended to the list, provided its £ 2 -distance to all orbits 
already on the list exceeded the threshold r = 10 -5 . 

There is a certain degree of arbitrariness in our choices of parameters. Regarding the 
choice of r and 5, recall that we currently have no means of predicting how close a (5-near 
biunimodular vector is to an actual biunimodular vector (if at all). Also, an a priori choice of 
the threshold r would have to be based on the knowledge of the local maxima of the mapping 
v i—^ ||Au||i on the torus. The above choice of r was motivated by the observation that the 
orbits of known biunimodular vectors in dimensions up to 13 had a minimal f 2 -distance 
> 10 -2 and 10~ 5 is safely below this value. 

The procedure to compute near biunimodular vectors was applied to all square-free num¬ 
bers between 2 and 15. In the following presentation of numerical results, we say that an 
orbit found by our algorithm is close to one from the literature, if the minimal distance 
(measured in || • || 2 ) between them is < 10 -5 . We observed that each time an orbit found 
by alternating projections was close to an orbit from the literature, the orbit cardinalities 
matched as well. 

• n = 2: The algorithm precisely found the one orbit consisting of the two vectors (1, ±i). 

• n = 3: The algorithm found one orbit of length 6, which was close to the orbit of the 
gaussian vector (1, w 2 ,1), where w = exp(27r?'/3). 

• 71 = 5: The algorithm found two orbits, each of length 10. Each orbit was close to one 
of the orbits represented by the gaussian vectors (1, w, w 4 , w 4 , oj) and (1, w 2 , w 3 , w 3 , w 2 ), 
where oj = exp(27rj/5). 

• n = 6: The algorithm found two orbits, one of length 12 and one of length 36. The 
orbit of length 12 was close to the gaussian vector, the orbit of length 36 was close to 
the vector described in ]M] formula (3.20)]. 

• n = 7: The algorithm found three orbits, with lengths 42,196 and 294, yielding 532 
solutions in total. The orbit of length 42 is close to the gaussian solution, whereas the 
orbit of length 196 is close to the Bjorck sequence (1,1,1, e* e , 1, e x6 , e ie ) mentioned in 
the introduction. The orbit of length 294 was close to the solution given in [34] formula 
(3.24)]. 

• n = 10: Here we obtained 10 orbits, with orbit lengths 20 (2 orbits), 200 (3 orbits), 
400 (4 orbits), and 800 (1 orbit), resulting in a total of 3040 solutions. For n = 10, the 
literature does not provide a complete list of solutions. 

• n = 11: Here 12 orbits were found, with lengths 110 (1), 484 (1), 1210 (7), 2420 (2) and 
4840 (1), resulting in a total number of 18744 vectors. Among these were the gaussians, 
all contained in the orbit of length 110. For n = 11, the literature does not provide a 
complete list of solutions. 

• n = 13: The algorithm found 20 orbits, with lengths 78 (2), 338 (2), 1014 (2), 1352 
(2), 2028 (7), and 4056 (5), yielding 40040 solutions in total. Gabidulin and Shorin [30] 
exhibited 25 orbits for n = 13, with 53222 solutions overall. Each of the orbits found 
via alternating projections was close to precisely one orbit representative from the list 
presented by Gabidulin and Shorin. The following orbits given in [30] were not found 
by the algorithm: An orbit with 1014 elements, corresponding to the case e = 6 in [30] 
Subsection 7.4], two orbits with 2028 solutions, corresponding to lines one and seven 
in the table in [30] Subsection 7.4], and the first two orbits listed for the case e = 12, 
each having 4056 elements. 

• 7i = 14: The algorithm found 39 orbits, with lengths 84 (1), 196 (2), 392 (1), 588 (3), 
784 (2), 1176 (12), 2352 (13), 4704 (5), resulting in a total of 72408 vectors. Among 
these were the gaussians, all contained in the orbit of length 84. For n = 14, the 
literature does not provide a complete list of solutions. 
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• n = 15: The algorithm produced 46 orbits, with lengths 900 (2), 1800 (21), 3600 (20), 
7200 (3), resulting in a total of 133200 vectors. The gaussian vectors are contained 
in two orbits of length 60, which were not found by the algorithm. For n = 15, the 
literature does not provide a complete list of solutions. 


Appendix A 


As we have already explained, Theorem 8.1 allows for an easy construction of all unitary 
matrices in the dyadic case U(2 n ). Recall, that every such matrix can be obtained from four 
matrices in U(2 n ~ 1 ). For n = 1 the procedure yields U(2) via (10). Here, we show how the 


procedure works in the case n = 2, to get a closed formula for 17(4). 


For a collection of 16 unimodular parameters a 3 , bj , Cj, z 3 £ T with j = 0,1, 2,3 we define 


the corresponding four matrices in U(2) via (10): 




ai(l + ao) 0103(1 — ao) 
a 2 (l-a 0 ) 0203(1 + a 0 ) 

Ci(l + C 0 ) CiC 3 (l-C 0 ) 
C 2 (l-c 0 ) c 2 c 3 (l + c 0 ) 




61 (1 + b 0 ) hb 3 (l-b 0 ) 
62(1 — 60) 6263(1 + 60) 

2 i(l + 2 0 ) 2 i 2 3 (1 - 2 0 ) 

-^2 (1 — ~Cl) ^2^3(l + “0) 


( 21 ) 


Theorem 8.1 asserts that every unitary matrix U £ 17(4) with entries ( U)k,i = Uk,i, where 


k,l = 1,2, 3,4, can be obtained as 


1 \A(I + Z) 

2 [B{I-Z) 


A{I - Z)C' 
B(I + Z)C 


This gives the following description of U (4): 


O 11 = h( a a i a 3 (1 — <* 0 ) 22(1 — z o) + s a i(l + ao)(l + §Al(1 + +)))) 

U21 = |(i 020 3 (l + ao) 22 (l - 2 o) + j 02 (l - a 0 )(l + § 2 l(l + 2 0 ))) 

O3I = i( — 36163(1 — 60)22(1 — 2o) + 361(1 + 6 0 )(1 — ^2i(l + 2 0 ))) 

041 = |( —36263(1 + 60)22(1 - 2 0 ) + 362(1 - 6 0 )(1 - 32i(1 + 2 0 ))) 

U12 = 3(301(1 + O 0 )2i 2 3 (1 - 2 0 ) + |Oia 3 (l - O 0 )(l + i2 2 2 3 (l + 2o))) 

022 = 3 ( 302(1 - O 0 ) 2 i 2 3 (l - 2 0 ) + |a 2 a 3 ( 1 + a 0 )(l + | 2 2 2 3 (1 + 2 0 ))) 

O 32 = 3( — 361(1 + 60)2123(1 - 2 0 ) + 36163(1 - &o)(l - 3 2 2 2 3 (1 + 2 0 ))) 

O 42 = h(~ 3 ^ 2(1 ~ 60 ) 2123(1 — 2 0 ) + 36263(1 + 6 0 )(1 — 3 2 2 2 3 (1 + 2 0 ))) 

O13 = 3 ( 301(1 + C 0 )(-iOia 3 (l - 00)22(1 - 2 0 ) + 303(1 + a 0 )(l - | 2 i(l + 2 0 ))) + 

302(1 - Co)(-iOi(l + 00)2123(1 - 2 0 ) + 30103(1 - a 0 )(l - 32 2 2 3 (1 + 2 0 )))) 
023 = 3(301(1 + Co)(-i02a 3 (l + a 0 ) 2 2 (l - 2 0 ) + 302(1 - a 0 )(l - |2l(l + 2 0 ))) + 
302(1 - c 0 )(—j0 2 (l - 00)2123(1 - 20) + ia 2 a 3 (l + a 0 )(l - |2 2 2 3 (1 + 2 0 )))) 
O33 = 3(301(1 + C 0 )(i6i6 3 (l - 60)22(1 - 2 0 ) + 1 61(1 + 60)(1 + |2i(l + 2 0 ))) + 
302(1 - C 0 )(j6i(l + 60)2123(1 - 2 0 ) + 36163(1 - 6 0 )(1 + |2 2 2 3 (1 + 2 0 )))) 

043 = 3 ( 301(1 + Co)(j&2& 3 (1 + 60)22(1 - 2 0 ) + 162(1 - 6 0 )(1 + i2i(l + 2 0 ))) + 
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jC 2 (l — C 0 )( j6 2 ( 1 — 60)2123(1 — ~o) + s6 2 6 3 (1 + 6 0 )(1 + i2 2 2 3 (l + 2g)))) 


M i4 — 5(^103(1 — co)(—jOia 3 (l — ao)2 2 (l — 2o) + Ja i(l + ao)(l — 323(1 + 20)))+ 
|c 2 c 3 (l + c 0 )(-|ai(l + a 0 )2i2 3 (l - 2 0 ) + iaia 3 (l - a 0 )(l - 3-223(1 + 2 0 )))) 
U24 = 3(3^1 ^3(1 — co)(— Ja 2 a 3 (l + ao)2 2 (l — 2 o) + |a 2 (l — ao)(l — |2i(l + 2o)))+ 
ic 2 c 3 (l + c 0 )(-ia 2 (l - 00)2323(1 - 2 0 ) + ia 2 a 3 (l + a 0 )(l - i2 2 2 3 (l + 2 0 )))) 
034 = |(|C!C 3 (1 - c 0 )(i6363(1 - 6 0 ) 2 2 (1 - 2 0 ) + 361(1 + 6 0 )(1 + |2i(l + 2 0 ))) + 
ic 2 c 3 (l + C 0 )( j6i(l + 60)2123(1 - 2o) + 56163(1 - 6 0 )(1 + | 2 2 2 3 (1 + 2 0 )))) 

044 = |(§cic 3 (l - c 0 )(j 6 2 6 3 (1 + 6o)2 2 (1 - 2o) + |6 2 (1 - 6 0 )(1 + §21 (1 + 2 0 )))+ 
3C2C 3 (1 + c 0 )( j 6 2 (1 - 60)2123(1 - 2 0 ) + 1 62 6 3 (1 + 6 0 )(1 + i2 2 2 3 (l + 2 0 )))) 


with 16 parameters dj, bj,Cj , Zj G T, where j = 0,1, 2, 3. 

By setting u±i = 1 in the above formula for U( 4) we will Hnd a description of 61(3), that 
is used in Theorem 6.3 to prove the existence of a biunimodular vector for every matrix in 
6/(3). This shall be executed as a series of observations. 

First, we set tin = 1 in the above formula for 6/(4). 


Observation A.l. tin := |Qaia 3 (l - a 0 )2 2 (l - z 0 ) + ioi(l + a 0 )(l+ i2i(l + 2 0 ))) = 1 iff 
a 0 = 01 = 2 0 = 21 = 1 . 


Proof. Clearly, tin = \{R,e + K), where e = (1,0), R = |(oi(l + ao), oia 3 (l — oo)) 
denotes the the first row of A from (21) and K = |( 2 i(l + z 0 ), 2 2 (1 — 2 0 )) denotes the first 
collumn of Z from therein. Thus, ||/ij ]2 = ||if|| 2 = 1. 

Since tin is an entry of a unitary matrix, we know that |un| < 1. This inequality can be 
proved directly in the following way 


l«n| = i \{R,e + K)\< \ (|(i?,e)| + |(i?,^)|) < \ (||£||||e|| + \\R\\\\K\\) = 1. 


Therefore, |un| = 1 means that \(R, e)| = ||i?||||e||, \{R,K)\ = ||f?||||A|| (so R = Ae, K = fie 
with A, /x G T) and | (R, e + K) | = |(i?, e)| + \(R, K)\ (so n = 1). Finally, itn = e + K) = 
1 gives A = 1. Hence, tin = 1 iff R = K = e, that is, do = 01 = 20 = 21 = 1. □ 

Then, we set ao = = Zq = Zi = 1 in the formula for 6/(4), to get a corresponding 

formula for 6/(3) given by the simplified entries Uf.^ with k : l = 2,3,4. By pulling out some 
factors from these entries we obtain the following matrix: //i 6 /( aj b jCj 2 )-D 2 , where D\ is the 
diagonal matrix D\ = diag(a 2 a 3 , 6363 , 6 2 6 3 ), _D 2 = diag(l, c 2 , c 3 ) and 


U, 


( a,b,c,z ) 


o±g) 

2 


(1 — Q )(l —c) 

4 


(l-q)(l+c) 

4 


(l-q)(l-b) (l+b)(l+c) (l+q)(l-b)(l-c) „ (l+b)(l-c) (l+q)(l-b)(l+c) 

4-44- 8 -4' 8 

(l-q)(l+b) „ (l-6)(l+c) (l+a)(l+b)(l-c) „ (l-iQ(l-c) (l+a)(l+b)(l+c) 

4 - 4”r 8 - 44- 8 


with a = 2 2 2 3 , 6 = 64 , c = C 4 and 2 = CiC 2 6 3 . This proves the following description of 6/(3) 

Observation A.2. 


6/(3) = {diag(\i, A 2 , ^z)U( a ,b,c,z) diag{l 1 A 4 , A 5 ) : a,b,c,z,Xj G T ,j = 1,...,5}. 
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By writing a = e l2a , b = e l2/3 , c = 
D[ = diag(e i “,e i (“ +/3 ), -ie^ a+ ^), D' 2 


e * 27 we get that U(a,b,c,z) = D'ffT^p^^D^, where 
= diag(l, e* 7 , —ie 11 ) and 


sin a cos 7 

z cos f3 sin 7 + cos a sin /? cos 7 , 
z sin /? sin 7 — cos a cos /? cos 7 

providing a parallel description of U (3): 

Observation A.3. 


-^(a>/3,7, 2 ) 


cos a — sin a sin 7 

— sin a sin /3 2 cos /? cos 7 — cos a sin ft sin 7 
sin a cos /? 2 sin /3 cos 7 + cos a cos ft sin 7 


17(3) = {diag(\ 1 ,\ 2 ,\3)T( ol ,p niZ )diag{l,\ i ,\ 5 ) : a ,£,7 £ [0,27t],z, Aj £ T,j = 1, ... ,5}. 


Appendix B 


In order to conduct a more detailed study of the alternating projection algorithm we begin 
by proving the following 

Proposition B.l. Let A £ U(n), v £ T n and Pa{v) = sign(A*(sign(Av))). Then, the 
condition ||Aw||i > n — 5 with 0 < S < ^ implies that 

(a) || 1 — |Au||| 2 < V25, 

(b) \Av\ > | and \A*(sign(Av))\ > 

(c) ||APa^Hi = ||Ai;||i iff Pav = v, 

(d) H-Pau — w||oo < 2y/n(\\ Ap Av\\i - ||A?;|| 1 )3. 

Proof. 

To see that \Av\ is close to 1 = (1,1,..., 1) we use ||v ||2 = \pn and observe that 

II 1 ~ 1Au111 2 = n - 2\\Av\h + \\Av\\l = 2(n - ||Au||i) < 25, 

proving (a). This already assures that \Av\ > \ for 5 < | (otherwise \ < ||1 — |Au|||| < 25, 
forcing 5 > |). 

Since \Av\ > 0, we easily see that ||sign(Ai;) — Au|| 2 = ||1 — |Aw||| 2 < y/25. Thus, 


|| A*(sign(Aa)) - v\\ 2 < V 26 . 

Therefore, by the same token as before, <5 < | implies that |A*(sign(Az>))| > 


To prove (c) we note that, by (15) and ( |16| ), 

||Au||i = (w,A*(sign(Av))) < ||A*(sign(Ai;))||i < ||A(sign(A*(sign(Az;))))||i = ||AP a u||i 


and, therefore, 


0 < ||A*(sign(Au))||i - (v, A*(sign(Av))) < \\APavWx - ||A?;|| 


( 22 ) 


Thus, ||APai;||i = ||At>||i implies that ||A*(sign(Ai;))||i = (v, A*(sign(At;)), what entails 
v = sign(A*(sign(At>))) = Pav, since v £ T n and |A*(sign(Ai;))| > 0. 

To show (d) we remark, that |A*(sign(Av))| > \ together with ||A*(sign(Au))||i < n and 
( [22| ) allows to apply Lemma B.4 below, proving (d). □ 

To show Lemma B.4 we start with an elementary observation. 

Lemma B.2. Let a, b £ C and 0 < e < 1. If a + b = 1 — e and |a| +16| = 1 then ||a| — a\ < ffe. 
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Proof. Since |1 — e — a| = 1 — |a|, we get that (1 — e)2Re(a) = (1 — e ) 2 — 1 + 2|a|, so 

||a| — a | 2 = 2|a|(|a| — Re(a)) = e|a|(2 — e — 2Re(a)) = e|a| ^1 + — 

what is clear for |a| > 1 , and follows from e < 2 |a| in the case when |a| < |. □ 

The next lemma can be treated as a solution to a certain random walk problem on the 
complex plane. 

Lemma B.3. Let u £ C n . If (it, 1) > 0, then 

IIM-tC, <IMIi(INIi -<«,!))■ 


-2M 

1 - e 


< e, 


(23) 


Proof. Let u = (ui,uz,... ,u n ) £ C“\ {0}. Due to rescaling we can assume that 
||it||i = 1. Moreover, the problem is insensitive to permuting the entries of u , so it is enough 
to show that ||iti| — u \\ 2 < 1 — (u, 1 ). The only hard part is to realize, that this estimate 
boils down to the previous lemma. 

Let us denote 1 — (u, 1 ) by e. Since 0 < (u, 1 ) < ||w||i = 1, it is clear that 0 < e < 1. 
Moreover, let us set a := u\ and x := Efc= 2 Ufc - Clearly, a + x = 1 — e. This means that 
the circle with the center at a and radius |.t| intersects the interval [ 0 , 1 ] at the point 1 — e. 


Since 


<EL 2 KI = i- 


the circle with the center at a and radius 1 — lal intersects 


the interval [0, 1 ] at some point 1 — e' with 0 < e! < e (this holds despite the fact, that the 
intersection can have two points). This proves an existence of b £ C such that a + b = 1 — e' 
and |a| + | 6 | = 1 with 0 < e 1 < 1, so Lemma B.2 yields ||a| — a | 2 < e' < e. □ 

Finally, we prove 

Lemma B.4. Let w £ C n and v £ T n . If 0 < - < \w\ and (v,w) > 0, then 


II sign{w) - v\\ OO — <VlMli(IMIi - (v,w)). 


(24) 


Proof. Clearly, 


:||sign(zo) - v\\oo < ||ra(sign(M;) - = |||«w| - vw\\ oo- 


Since for u := vw we have (u, 1) = ( v,w ) > 0 and 


Lemma B.3 yield (241 immediately. 


Mil = IMIli, the above estimate and 

□ 


After proving Proposition |B. l| we can take on understanding the workings of the algorithm 
searching for the biunimodular vectors of a given A £ U(n). Recall that we start with an 
arbitrary vector Vq £ T" and successively apply Pa (given in Proposition B.l), so that 
Vj = P j a Vq, until we reach a J such that 


||AVj||i >n-8>n-\ 

with a certain precision 0 < S < i. As we mentioned in Subsection 


(25) 

there is no guarantee 

that for a given starting vector Vq we will find J £ N so that (|25| holds. Hence, in practice, 


10.1 


the algorithm has to be applied to many random starting vectors in T". 

Thus, to facilitate the further discussion of the algorithm, we will concentrate on the 
starting vectors Vq satisfying the initial condition 


Will > n- 


(26) 


Under this assumption we can gather our findings in the following 

Proposition B.5. Let A £ U(n), Vq £ T n and Vj = P a Vq for j £ N. If HAVoH! > n — 
then for all j £ N U {0} we have that 

(a) \\AVj\h < ||AV J+1 || 1 < n, 
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1 and \ A*(sign(AVj))\ > 1 


(b) V 3 - G T n , \AVj\ > 2 

(c) = \\AVjh iffV j+1 = V j , 

(d) \\Vj + i - VJIloo < 2 y /n(\\AV j+1 \\ 1 - ||^ 

(e) lim 


2 ’ 


l j-HX || Vj +1 


Vi+l - Vi II no =0. 


Proof, (a) recalls ( [T7| ) a nd | IA|| 


oo—> 1 

Vj G T" then, by Proposition B.l'b), \AVj\ > 


the latter, that is Vj+\, must belong to T 
(e) follows from (d), since (a) assures that the sequence (||AV) 


(b) is a matter of a simple induction. If 

so the signum of 
B.ll And 


n. 

§ and |A*(sign(AP,-))| > \ 

(c) and (d) follow from Proposition 


is convergent. 


□ 

Since the sequence (||AV)j|i)j e N is convergent one can hope, that Vj shall converge to a 
fix point of Pa- This, however, is hard to assure even under the initial condition ( |26[ ). This 
lack of proof for converge nce can be traced to the square root appearing in Lemma ]B.2| that 
is hard to remove, since (23) gives ||a| — a| > \Je\a\ for |a| < Therefore, we can only fall 

Vj+1 


on linij^oo 


-v,\ 


= 0 to draw a standard conclusion 


Corollary B.6. If the initial condition (26) holds, then the cluster points of (Vj)j^ belong 
to the set Sa '■= {u G T" : Pav = v}. Moreover, if S a is finite, then Vj converges to a point 
in Sa- 

Proof. Here, we use the notation provided in the two previous propositions and by 
convergence we mean the convergence in any F'-norm, since they are all equivalent. 

Let L denote the limit of If a subsequence (Vj k )ke n of (Vj)j e ^ converges to 

V, then PAVj k converges to PaV, since by (b) of Proposition 


B.5 


the sequence (Vj)jgN is 
contained in a closed subset of an open domain, where Pa is continuous. Thus, ||AP||i = 
lim^oo ||AVjf fc ||i = L and ||APaV||i = lim^^oo \\AP A V jk \\i = lim^oo ||AV'^ + i||i = L, so 


||AF||i = HAP^PIIi. Since V G T" and ||AP||i > n — g, Proposition B.l'c) implies that 
P A V = V, hence V G Sa- 

If Sa is finite, then its points are separated by open balls. Thus, almost all points of 


are contained in the union of these balls. By Proposition B.5 e), there is only one 


of these balls containing almost all points of (V^eN, what assures the convergence of Vj. 


□ 


Since Bi(A) C Sa, the finiteness of Sa is very hard to establish and, in the Fourier case, 
brings us back to Conjecture o 
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